Preprint typeset in JHEP style - HYPER VERSION 



TAUP - 2825-06 

February 2, 2008 
hep-ph/0606260 



The BFKL Pomeron Calculus: Probabilistic 
o Interpretation and High Energy Amplitude. 

bJO 

^ M. Kozlov *, E. Levin^ and A. Prygarin* 

Department of Particle Physics, School of Physics and Astronomy 
^NJ Raymond and Beverly Sackler Faculty of Exact Science 
Tel Aviv University, Tel Aviv, 69978, Israel 

\o '■ 
(N ■ 

\^ Abstract: In this paper we continue to pursue a goal of finding the effective theory for high energy in- 
teraction in QCD based on the colour dipole approach, for which the BFKL Pomeron Calculus gives a low 
^) energy limit. The two key problems, that we consider are the following: the probabilistic interpretation 
of the BFKL Pomeron Calculus and the possible scenario for the asymptotic behaviour of the scattering 
^p-iamplitude at high energy in QCD. In this paper we show that the generating functional approach is equiv- 
Qnalent to the BFKL Pomeron Calculus and leads to a clear interpretation of this calculus as an alternative 
r"| iJescription of the system of interacting colourless dipoles. 

j>! We calculated the two BFKL Pomerons into one BFKL Pomeron vertex directly in the dipole approach and 
show that this merging can be described as the decay of two dipoles into four dipoles in the dipole approach. 

5_( This result means the JIMWLK approach can give us all needed vertices for the BFKL Pomeron calculus 
■ — and, therefore, they together give the selfconsistent and complete theoretical approach to high energy 
scattering in QCD. 

In this paper, we developed the semi-classical approach to the functional evolution equation for the gener- 
ating functional. Using this method we found the solution in the entire kinematic region starting from low 
energies for the simplified case of the toy model in which we assume that interaction of dipoles does not 
depend on their sizes. In the general case we find the asymptotic solution at ultra high energies as well as 
the first correction to this solution. 

We considered the semi-classical approach to the BFKL Pomeron Calculus and recovered the Mueller- 
Shoshi band for the kinematic region where we can trust the linear BFKL equations. The solution for the 
scattering amplitude deeply in the saturation region was found and discussed. 
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1. Introduction 

The simplest approach that we can propose for high energy interaction is based [|^, ^] on the BFKL Pomeron 
and reggeon-like diagram technique for the BFKL Pomeron interactions @, |5|, 0. This technique , 
which is a generalization of Gribov Reggeon Calculus |^], can be written in the elegant form of the functional 
integral (see Ref. and the next section); and it is a challenge to solve this theory in QCD finding the high 
energy asymptotic behaviour. However, even this simple approach has not been solved during three decades 
of attempts by the high energy community. This failure stimulates a search for deeper understanding of 
physics which is behind the BFKL Pomeron Calculus. On the other hand, it has been known for three 
decades that Gribov Reggeon Calculus has intrinsic difficulties Q that are related to the overlapping of 
Pomerons. Indeed, due to this overlapping we have no hope that the Gribov Reggeon Calculus could be 
correct in describing the ultra high energy asymptotic behaviour of the amplitude. The way out of these 
difficulties we see in searching for a new approach which will coincide with the BFKL Pomeron Calculus 
at high, but not very high, energies (our correspondence principle) but it will be different in the region of 
ultra high energies. In a spirit of the parton approach we believe that this effective theory should be based 
on the interaction of 'wee' partons. We consider, as an important step in this direction, the observation 
that has been made at the end of the Reggeon era 0, |ll|, H], that the Reg geon Calculus can be reduced 



to a Markov process for the probability of finding a given number of Pomerons at fixed rapidity Y. 
Such an interpretation, if it would be reasonable in QCD, can be useful, since it allows us to use powerful 
methods of statistical physics in our search for the solution. 

The logic and scheme of our approach looks as follows. The first step is the Leading Log (1/x) 
Approximation (LLA) of perturbative QCD in which we sum all contributions of the order of {as ln(l/x))"'. 
In the LLA we consider such high energies that 

as ln{l/x) w 1; while as <C 1 (1.1) 

It is well known that the LLA approach generates the BFKL Pomeron (see Ref. Q and the next section) 
which leads to the power-like increase of the scattering amplitude {A oc ^^^„Jo,,^^o) with 
u){n = 0, = 0) oc as)- 

The second step is the BFKL Pomeron Calculus in which we sum all contributions of the order of 

(a^ 1 \^ o? 1 

l3y — 1 — 1\ — ml therefore — — - — -r- 1 (1.2) 

where as = osNc/tt. 



The structure of this approach as weh as its parameter has been understood before QCD [14| and was 
confirmed in QCD (see Refs. fl], |, |,|, |, 0, |l6|). This calculus extends the region of energies from 
ln(l/x) ~ l/as of LLA to ln(l/j;) ~ (l/as) ln{N^ /as)- The BFKL Pomeron Calculus describes correctly 
the scattering process in the region of energy: 

1 fN^\ 11 , , 

— In ^ ^ In- ^ 1-3 

as \ as J X 

For higher energies the corrections of the order of (q;| ln(l/x))" should be taken into account making all 
calculations very complicated. 

Our credo is that we will be able to describe the high energy processes outside of the region of Eq. (|1.3|), 
if we could find an effective theory which describes the BFKL Pomeron calculus in the kinematic region 



given by Eq. (O), but based on the microscopic degrees of freedom and not on the BFKL Pomeron. In 
so doing, we hope that we can avoid all intrinsic difficulties of the BFKL Pomeron calculus and build an 
approximation that will be in an agreement with all general theorems like the Froissart bounds and so on. 
Solving this theory, we can create a basis for moving forward considering all corrections to this theory due 
to higher orders in as contributions, running QCD coupling and so on. 

The goal of this paper is to consider two main problems: the probabilistic interpretation of the BFKL 
Pomeron Calculus based on the idea that colour dipoles are the correct degrees of freedom at high energy 
QCD ||l^; and the possible solution for scattering amplitude at high energy. We believe that colourless 
dipoles and their interaction will lead to a future theory at high energies which will have the BFKL Pomeron 
Calculus as the low energy limit (see Eq. ( |1.3D ) and which will allow us to avoid all difficulties of dealing 
with BFKL Pomerons at ultra high energies. 

Colourless dipoles play two different roles in our approach. First, they are partons ('wee' partons) 
for the BFKL Pomeron. This role is not related to the large Nc approximation and, in principle, we can 
always speak about probability to find a definite number of dipoles instead of defining the probability to 
find a number of the BFKL Pomerons. The second role of the colour dipoles is that at high energies we can 
interpret the vertices of Pomeron merging and splitting in terms of probability for two dipoles annihilate 



in one dipole and of probability for decay of one dipole into two ones. It was shown in Ref. [17| that 
P 2P splitting can be described as the process of the dipole decay into two dipoles. However, the 
relation between the Pomeron merging (2P P) and the process of annihilation of two dipoles into one 
dipole is not so obvious and it will be discussed here. 

This paper is a next step in our programme of searching the simplest but correct approach to high 
energy scattering in QCD in which we continue the line of thinking presented in Refs. ||l^, ^ 21 1. The 
outline of the paper looks as follows. 

In the next section we will discuss the BFKL Pomeron Calculus in the elegant form of the functional 
integral, suggested by M. Braun about five years ago [H. We will show that the intensive recent work on 



this subject |22, |2^, confirms the BFKL Pomeron Calculus in spite of the fact that these attempts were 
based on slightly different but not more general assumptions. 

In the third section the general approach based on the generating functional will be discussed. The set 
of equations for the amplitude of n-dipole interaction with the target will be obtained and the interrelations 



between vertices of the Pomeron interactions and the microscopic dipole processes will be considered. It 
will be shown that the generating functional approach is equivalent to the BFKL Pomeron Calculus in 
the kinematic region of Eq. ( |1.3| ) and leads to a clear interpretation of the BFKL Pomeron Calculus as an 
alternative description of the system of interacting colorless dipoles. 

The Fourth section is devoted to the probabilistic interpretation of the BFKL Pomeron Calculus and, 
in particular, the process of merging of two Pomerons in one. We calculate the vertex of this merging 
directly in the dipole approach and show that it is related to the process of two dipole decay into four 
dipoles. We discuss the toy model which simplify the QCD interaction and allows us to see the main 
properties of high energy amplitude. In particular, we are going to discuss the solution of the equations 
for the asymptotic behaviour of the scattering amplitude at high energies. We find the solution for the 
generating functional in the entire phase space of both variables: rapidity Y and the variable u which is a 
conjugate variable to the number of dipoles. 

In the fifth section we discuss the solution to the generating functional at high energies. We confirm 



the analysis of Ref. |21] and find the asymptotic solution suggested there. However, we study in more 
details the mechanism of approaching the ultra high energy limit for the scattering amplitude. 

In conclusion we are going to compare our approach with other approaches on the market. 
2. The BFKL Pomeron Calculus 

2.1 The general structure of the BFKL Pomeron calculus. 

We start with a general structure of the BFKL Pomeron calculus in QCD. The BFKL Pomeron exchange 
can be written in the form (see Fig. |T|-1) 



_2 

^(Fig.|-1) = K(g)Gp(ri,r2;6|yi-l2)(g)Ki « ^ exp (u;(n = 0, = 0) - Ya)) (2.1) 

with Lo{n = 0, = 0) oc as and denotes the all needed integrations. — I2 = ln(l/x). To understand 
the main parameters of the BFKL Pomeron calculus, it is enough to compare the contribution of the 
first 'fan' diagram of Fig. 1-2 with the one BFKL Pomeron exchange. This diagram has the following 
contribution 



A(Fig. ^2) = r dY' Vu (g) Gp(ri, /; b\Yi - Y') (g) r(l ^ 2)G|,({r'}, rs; h\Y' - Y^) (g) (2.2) 

JY2 



oc ^" ^ ^1 exp (2 a;(n = 0, 1/ = 0) Y) oc exp (2 c^(n = 0, = 0) Y) 

u}(n = 0, = 0) 

where Y = Yi — Y2] Ri and r2 are the sizes of the projectile and target dipoles while {r'} denotes all dipole 
variables in Pomeron splitting and/or merging. 

_2 

One can see that the ratio of this two diagrams is proportional to -r^ exp (a;(n = 0, = ^){Yi — I2)) 
which is the parameter given by Eq. (1^2). When this ration is about 1 we need to calculate all diagrams 



A 



with the Pomeron exchange and their interactions (see Fig. Ill-a - Fig. |- f ). Ah vertices, that are shown in 
Fig. |l[ has been calculated in Refs. |Q, |5| and they have the following order in 05 ^: 

-2 -2 

u;(n = 0,J^ = o) oc as; r(l ^ 2) oc ^ ; r(2 ^ 1) oc ^ ; 

r(2^2) « ^;r(2-3) oc ^;Koc ^;^dOC (2.3) 
Using Eq. (^|^) we can easily estimate the contributions of all diagrams in Fig. |^. Namely, 



A(Fig. 0-3) oc Vu Vd ^^l^ ^^^^^^ J.^ exp (2a;(n = 0, 1/ = 0)(yi - Y^)) 

(-2 \ 2 

^ exp(u;(n = 0,z. = 0)(yi-y2))J = L\Y) ■ (2.4) 

^(Fig. 0-4) oc K t^^^^'^^ ^ 2) (Fi - Y^) exp (2u;(n = 0, 1/ = 0)(yi - Y^)) 

oc ^(Fi-y^) (^exp(a;(n = 0,.. = 0)(yi-y2))) =^yL2(y); (2.5) 
^(Fig. 0-5) oc Vu V! ^^Ir^^ll^^^Zf) = 0' = - ^2))) 

1 / ~^ \ ^ 1 

A(Fig. 0-6) oc Vu Vd 2) r(2 ^ 1) ^ ^ _ ^ = 0, 1/ = o)(yi - Y2))) 

LV^yn = U, 1/ = (Jj 

oc ^(Fi-y^) (^exp(a;(n = 0,r. = 0)(yi-y2))) =^YL\Y)- (2.7) 

with L(y) = (a|/Ar2) exp (w(n = 0, = 0) y). 

As we have mentioned in the introduction the BFKL calculus sums all diagrams at so high energy that 
parameter L{Y) is of the order of 1 (see Eq. ( |l.2D ). In this kinematic region we need to take into account 
the diagrams of Fig. -1 , Fig. |ll-2 and Fig. |l]-3 (see Eq. (^rT|),Eq. (|2l^ and Eq. (^^). Indeed, diagrams 
of Fig. [11-4 and Fig. 0-6 ( see Eq. (U) and Eq. (|3)) are small since {al/N"^) Y <^ l/N"^ «; 1 in the 
kinematic region of Eq. (|1.3| ), while the diagrams of Fig. 0-5 (see Eq. ( |2.6D ) is small at L (Y) « 1. 

The first conclusion that we can derive from this analysis that in the kinematic region where L {Y) ~ 1 
we need to take into account all diagrams with r(l 2) and r(2 — > 1) vertices while the diagrams with 
r(2 2) and r(2 3) vertices give small, negligible contributions. 

However, if L (y) oc Nc one can see from Eq. - Eq. (^ that all dia grams give so essential 

contributions that we have to take them into account. For what follows it is interesting to notice that the 
vertex r(2 ^ 2) oc a^/N^ can be neglected even at such large values of L (Y). 

^In Eq. ( ^.3| ) we use the normalizations of these vertices which are originated from calculation of the Feynman diagrams. 
In the dipole approach we use a different normalization (see below section 3 and 4) but all conclusions do not depend on the 
normalization. 



r 



Vu Vd r(i->2) r(2->i) r(2->2) r(2->3) 




Figure 1: The BFKL Pomeron interactions and the examples of the diagrams of the BFKL Pomeron calculus in 
QCD. The solid line describes the Pomeron exchange while the double line stands for the dipole. 



Finally, we can conclude that the first step of our approach can be summing of the diagrams with 
r(l 2) and r(2 1) vertices. Namely, this approach we discuss below in details. 

However, we would like to stress that we need to make an additional assumption inherent for the 
BFKL Pomeron calculus: the multi-gluon states in t-channel of the scattering amplitude lead to smaller 
contribution at high energies than the exchange of the appropriate number of the BFKL Pomerons (see 
more in Ref . |2l|] ) . The real argument, supporting this assumption, stems from the numerous attempts to 
find the intercept of these states, larger than the one for the multi-Pomeron exchanges |25| that failed. 



2.2 The path integral formulation of the BFKL calculus 

The main ingredient of the BFKL Pomeron calculus is the Green function of the BFKL Pomeron describing 
the propagation of a pair of gluons from rapidity Y' and points x'l and X2 to rapidity Y and points xi 
and X2 ^ Since the Pomeron does not carry colour in t-channel we can treat initial and final coordinates as 
coordinates of quark and antiquark in a colourless dipole. This Green function is well known |2^ and has 
a form: 

G{xi,X2;Y\x[,X2;Y') = Q{Y -Y')x (2.8) 

00 „ 

X du d^xoe'^^''''''^^^-^'^ X{n,u)E{xi,X2;xo\u)E*{x[,x'2;xo\u) 

n=—oo 

^Coordinates xt here are two dimensional vectors and, strictly speaking, should be denoted by Xi or Xj. However, we will 
use notation Xi hoping that it will not cause difficulties in understanding. 




Figure 2: The graphic form of tlie triple Pomeron vertex in the coordinate representation. 



Eix,,xr,xoW) = (-^ (2.9) 



where vertices E are given by 



3^10 2^20/ \a^io^20 

where Xik = Xi — Xk, Xi = Xi,x + ixi^y ^ ,x* = Xi^x + ixi^y ; h = {I — n)/2 + iu and h = 1 — h*. The energy 
levels uj{n, v) are the BFKL eigen values 

u;(n, u) = as (^{l) - Rei; + ii^^ ) (2.10) 

where il^{z) = dlnT{z)/dz and T{z) is the Euler gamma function. Finally 

= [(^ + 1)2 + 4^,2] _ 1)2 + 4^2] (2-11) 

The interaction between Pomerons is depicted in Fig. ^ and described by the triple Pomeron vertex 
which can be written in the coordinate representation ^ for the following process: two gluons with 
coordinates x'l and X2 at rapidity Y"/ decay into two gluons pairs with coordinates and Xg at rapidity 
Y2 and X2 and X3 at rapidity due to the Pomeron splitting at rapidity Y. It looks as 



7ra| f (fxi(fx2(fx;>, . I I v\*r \ ro1o^ 

2^ / ^2 ^2 ^2 — (G(xi,X2;yi|xi,X2;y)tLi^2)- (2-12) 



;,2 /■ ^2^_^ ^2 ^2^ 

^2 ^2 ,y.2 

''^12 -^^23 '''13 

G(xi, X3; x'3'; ^2') G(X3, X2; x'2'; ^3') 

where 

T^i,2 = rf^pjpl with / = -V^ (2.13) 



^Xi^x and are components of the two dimensional vector Xi on x-axis and y- axis 



and arrow shows the direction of action of the operator L. For the inverse process of merging of two 
Pomerons into one we have 

-K a% f d?xi (P X2 (Pi 



Nr. 



I d xid X2 d X3 Qi^^//^^ii^. Y(\xi,x3; Y) G{x'^, 4'; ^^2^3, X2; Y) ■ 

■J ^12 ^23 ^13 

(ri2G(xi,X2;y|x;,x'2;y/)) (2.14) 

The theory with the interaction given by Eq. (|2l^ ) and Eq. (|2l|) can be written through the functional 
integral 

= j D^D<^>+e^ with S = So + Si + Se (2.15) 

where 5*0 describes free Pomerons, Si corresponds to their mutual interaction while Se relates to the 
interaction with the external sources (target and projectile). From Eq. ( 2.12| ) and Eq. (2.14) it is clear that 

So = j dY dY' d^xid^X2d^x\d^X2^^{xi,X2]Y)G~^{xuX2]Y\x^,X2;Y')(^{x^,X2-X) (2.16) 
Si = ^^ [ dY [ d'Y^'^'''' ■m,2Hxi,X2;Y)) . cI>+(xi,X3;y)$+(x3,X2;y) + h.c.} (2.17) 

i\c J J X^2 2^23 -^13 

For Se we have local interaction both in rapidity and in coordinates, namely, 

Se = - j dY d^Xid^X2{^{xi,X2]Y)Tpr{xi,X2]Y) + ^+{xi,X2\Y)Ttar{xi,X2]Y)] (2.18) 

where r^^ {Ttar) stands for the projectile and target, respectively. The form of functions r depend on the 
non-perturbative input in our problem and for the case of nucleus target they are written in Ref. Q. 

For the case of the projectile being a dipole that scatters off a nucleus the scattering amplitude has 
the form 

T{x„y,;Y) = \x^,y,;Y) = - j eSi^M]\,^^, ^^'''^ 

where extra as comes from our normalization and where we neglect term with Tpr in Eq. ( p. 18 ). 

Generally, for the amplitude of interaction of n dipoles at rapidity Y we can write the following 
4 



expression 



TH(x,.„...x.,..;y) = (-ir(^j ^ (2-20) 

The extra factor (—1)" stems from the fact that in Se the source for a projectile as well as for a target, 
has an extra minus sign. 

It is useful to introduce the Green function of the BFKL Pomeron that includes the Pomeron loops. 
This function has the form 

r^i ^1 / D^D^^^^{x^,yv,Y)^{x2.y2;Y') e^t^-'^+l 



''starting from this equation we use notations Xi for the coordinates of quark while j/i denote the coordinates of antiquarks. 
For rapidity we will use Y . 



For further presentation we need some properties of the BFKL Green function [26|: 
1. Generally, 



G-'{xi,X2;Y\x[,x'^-Y') = p\pI 1^-^ + = {J^ + plpl (2.22) 
Hf{x^,X2;Y) = ^ [ ^-^^ {f{xi,X2;Y) - f{xi,xs;Y) - f{x,,X2;Y)) ■ (2.23) 



'^X^ Xi. 
^^23 ^13 

2. The initial Green function (Gq) is equal to 



2 2 2 2 

Go{xi,X2;Y\x'i,X2;Y) = vr^ln^j^ In — (2.24) 

1,2' ■^l',2 ■''1,2 -^l',2' 

This form of Go has been discussed in Ref . p6[ . In appendix A we demonstrate that this expression for 
Gq{xi,X2] Y\x'i,X2; Y) stems from u = uj{n = 0, v) term in sum of Eq. (p.SD- Only this term is essential at 
high energies since all other terms lead to decreasing with energy contributions. 

3. It should be stressed that 

VlvlGo{xi,X2-^Y\x\,x'2;Y) = (2.25) 

(2 irf ((^(2) _ ^/^) ^(2) _ 4) + 5(2) {^xi - x'2) 5^'^'^ {X2 - x'l) ) 



4. In sum of Eq. ( p. 8]) only the term with n = is essential for high energy asymptotic behaviour since 
all u;(n, u) with n > 1 are negative and, therefore, lead to contributions that decrease with energy. Taking 
into account only the first term one can see that G is the eigen function of operator L13, namely 

Li2G(xi,X2;y|x;,4;y') = -J--G(xi,X2;r|x;,x'2;n « G(xi, X2; y|x;, x'2; F'); (2.26) 

The last equation holds only approximately in the region where v 1, but this is the most interesting 
region which is responsible for high energy asymptotic behaviour of the scattering amplitude. 

2.3 The chain of equations for the multi-dipole amplitudes. 



Using Eq. ( p. 15 ) and Eq. ( p.l9| ) we can easily obtain the chain equation for multi-dipole amplitude T^") 



noticing that every dipole interacts only with one Pomeron (see Eq. ( p.l9| )). 

These equations follow from the fact that a change of variables does not alter the value of functional 
integral of Eq. (|2.15| ). In particular, Z[<I>,$+] = <I>'+](see Eq. (|2.15| ) ) where <!>'+ = + e{x,y) 



with small function e{x,y). Therefore, 

j i?$i^$+e^[*'*^] = j Z)$i^$'+e^[*'*'+l (2.27) 
Substituting = + e(x,z/) and expanding this equation to first order in e, we find 



+ 



= J D^D^+e^l*'*""] X (2.28) 

J dY dY' d^xi d^X2 d^X2 e{xi,X2;Y) G-^{xi,X2;Y\x[,X2;Y') Xg; 1^') 
d'^xi d^X2 d^x-i 



2 VTQ 



dY 



111 

/y*iJ rr**^ /ya— ^ 

■^Vl -^23 -^13 



{(e(xi,x2;y)ri^2) • ^(xi,X3;y)^(x3,X2;y)+ 



2{Tx2^{x\,X2\Y)') ■ e(xi,X3;y) ^>+(x3,X2;l")} - j dY Sx\Sx2€{xx,X2\Y)Ttar{x\,X2\Y)\ 

We redefine the integration variables in the third term as follows 

2 {T^^2^[xx.X2-, Y) ) • e(xi, X3; Y) $+(x3, X2; ^ 2 (ri_3$(xi, X3; ) • e(xi, ^2; $+(x2, X3; 



Using the expression for the Hamiltonian Eq. 2.22 and the Casimir operator Eq. 2.13 we define a new 
variation parameter e{xi,X2',Y)p1p2- In terms of this variation parameter Eq. 2.25 reads as 



= / D^D^+e^l*'*""] 



(2.29) 



+ 



[j dYd^xid^X2 e{xi,X2;Y)plpl(^-^ + H^ <!>{xi,X2;Y) 

^■^"l f y^r f d^Xid^X2d^X-i , 2 2\ ^, . 

^ ' dY I — ^ — ^ — ^ — ■{{e{xi,x2;Y)pipi) ■ $(xi, X3; y) $(^3, X2; y)+ 



^2 ^2 0^2 
''^12 -^^23 -^13 



2 2 r 2 2 

2(ri3$(xi,X3;y)) • e(xi,X2;y)^^+(x3,x2;y)}- / dY d^xid^X2e{xi,X2;Y)^nar{xi,X2]Y)] 

P1P2 J P1P2 

We denote the new variation parameter by e{xi, X2;Y) = p1p2e{xi, X2;Y) and use the property of the 



initial Green function Eq. 12.251 to rewrite -^-^ in terms of Gq as follows 

' ' P1P2 



-^e{xi,X2;Y) 
P1P2 



(2vr)^ 



Go(xi, X2; y|x'^, y)e(x'^, Y)d'^x'id'^x'2 



Thus, Eq. 2.29 can be written as 



+ 



= J DcI>Z)cI>+e^[*'*''] X 
[j dYd^xid''x2 eixi,X2;Y)(^-^ + H^^xi,X2;Y) 
27ra| f f d^x^d^dy^ .^~^^^^^^^.Y) cD(xi, X3; F) a>(x3, X2; y)+ 



(2.30) 



(2.31) 



/y"2 ^2 ,y.2 

-^12 -^^23 '''13 



2(ri_3$(xi,X3;y)) • 1^ j Go{xi,X2;Y\x[,x'2;Yre{x[,x'2;Y)d^x[d''x'2^ <D+(x3, X2; y)} 

-/ 



dY d xid X2 



{ (2vr)^ 



Go{xi,X2;Y\x[,X2;Y)e{x[,x'2;Y)d^x[d'^X2 } Ttar{xi,X2;Y)] 



Noting that the r.h.s. of Eq. ( p.31|) should vanish for any possible variation e{xi,X2; Y) we obtain 

= j Z)$D$+e^[*'*^] x[(^-^ + H^'^ixi,X2;Y) (2.32) 



+- 



2Trai 



+ 
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(fx[ (fx'2 d^x-i , 4 



^f^$(xi, X3; Y) $(X3, X2; Y)+ 

^23 ^13 



,2 ^2 ^2 
2"' 

1 



ry*^ ry^^ ry^^ 



{xvi X3; Y) ) Go x'2; Y\xi,X2]Y) $+(x3, 4; y)} 



(2vr)^ 



d X^d X2GQ{Xi,X2;Y\xi,X2]Y)Ttar{Xi,X2:,Y)] 



where in the third and last terms we interchanged (x'^jXg) ^ (xi,X2). We notice that the third and last 
terms are independent of rapidity and can be absorbed in the initial condition. It is obvious for the last 
term which represents the target source. To show this for the third term we use the property of the Casimir 
operator at high energies (n = , = 0) 

ri 3$(xi,X3;y) ^ $(xi,X3;y) 



and the definition of the Green function(see Eq. 2.21 ). We see that the third term results into the product 
of two initial Green functions which are independent of rapidity. 



Now we can use the definition of the amplitude defined in Eq. |2.19| and Eq. p.20| to rewrite Eq. |2.32 
in a simple form 



arW(xi,x2;y) 



dY 



as 
2tt 



cf z K{xi, X2; z) 



(2.33) 



(r«(xi,z;y) + r«(z,X2;y) - r«(xi,X2;y) - t'-^\Y-xuz-z,X2;Y)) 
where kernel K{x,y\z) is defined as 



K{x,y\z) 



{x - yf 



{x — zY (z — yY^ 



(2.34) 



and the Hamiltonian is given by Eq. 2.25 . 

This equation has a very simple meaning that is clear from Fig. ^. 
Starting from equation 

J L»$Z:i$+$(y;x4,X5) e^[*'*^l = j L> $ L> ^>(y; X4, X5) e 

we obtain the equation for the amplitude T^'^\ namely, 

5r(2)(xi,x2;x3,X4;y) as , ,^ , . 

— a z K{xi,X2\z) 



(2.35) 



dY 2-K 

{xi,z;x3,X4;Y) + T^'^\z,X2;x3,X4^;Y) - T^'^\xi,X2;x3,X4;Y) 



(2.36) 
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Figure 3: The graphic form of equations for the multi dipole amphtude. 



T^^''{xi,z;z,X2;x3,X4;Y)] + 



where function r2^i is equal to 



+ 2a|— ^ / (f X x" r2-»i(xi, X2; X3, X4|x', x") V^/V^.//T^^^(x', x" \ Y)\ 



^2-^\{xx,y\\X2,yi\x,y) = / (fzK{x,y;z)Go{xi,yi\x,z)Go{xi,yi\z,y) 



(2.37) 



Deriving Eq. ( ^33] ) and Eq. (|23^ ) we use Eq. ( p. 24] ) and Eq. ( |2.25| ) as well as normalization condition (see 
Eq. ( 2.19D ) for the scattering amplitude. These two equations are the same as in Ref. [|7|- This shows that 
in the papers |22, B^, actually, the same approach as the BFKL Pomeron Calculus has been developed 



(much later, of course), in spite of the fact that authors think that they are doing something more general. 

Assuming that T^^) = T^^) T^^) we obtain the Balitsky-Kovchegov equation ||2^, |2^. We can do this 
only if we can argue why the Pomeron splitting is more important than the Pomeron merging. For example, 
this assumption is reasonable for scattering of the dipole with the nucleus target. Generally speaking, the 
splitting and merging have the same order in as ( see Eq. {^J^ and Eq. (|1|). In Eq. (p^) and 
Eq. ( 2.36| ) these two processes look like having a different order of magnitude in 05 but this fact does 
not interrelate with any physics and reflects only our normalization. However, we will see that for a 
probabilistic interpretation the correct normalization is very important. 



3. Generating functional and probabilistic interpretation 

3.1 Statistical physics analogy: Langevin equation and directed percolation. 

The functional of Eq. ( p. 15 ) is reminiscent of the partition function of statistical mechanics. Indeed, the 



1 n 



partition function has a general form 

Z[H] = e-^^[^l = ydse-w/rf-W^)-^^^W) (3.1) 



where F[H] is the Helmholtz free energy. As an example, Eq. ( |3.1| ) is written for the system of spins with 
the energy density 'H(s) in the external magnetic field H. The integration is over all possible configurations 
of spins in the system. 

Comparing Eq. (jsj) and Eq. one can see that Eq. ( |2.15| ) describes a statistical system with 

kT = 1 and with F[—Ttr\ = —S. The form of 5^ suggests that H = —Ttr plays a role of the external field. 

Eq. ( 2.15| ) on the lattice looks as follows 



= j \{D^i{xi,yi-Yi)D^t{xuyi-^Yi)e^^'^-'^t\ (3.2) 



We can linearize the action using the Stratonovich transformation ||30[] but first we simplify the ex- 
pression for action S introducing a new field {xi,X2]Y) = p\p\<^'^ {xi,X2]Y) and using Eq. ( 2.22 ) , 



(3.3) 



Eq. ( |2.26| ) and Eq. ( 2.30 ). In new notation the action has the form 

5o = y dY dY' d\<fyd\' (fy' {x,y;Y) {J^ + Hx',y';Y') 

Sj = (^J dYd^xd^yd^zK{x,y\z)<!>+'{x,y;Y)^{x,z;Y)^{z,y;Y) + (3.4) 

1 A ^2 ,2 ^2~>iixi,yi;,x2,y2\x3,y3) +', v\ff,+'f v\^/ v\\ 
j^-p^ [[d Xid yi _ ^^^4 [xi,yi;Y)^^ [x2,y2]Y)^{x3,y3]Y)\ 

while Se has the form 

5^; = -^ I dY d^xd^yd^x'd^y' ^+'{x,y;Y)Go{x,y\x',y')T{x',y';Y) (3.5) 

We demonstrate the idea of the statistical interpretation of our theory using the simplified version, so 
called toy model, in which we assume that there is no dependence of the interaction on the coordinate of 
dipoles (see Ref. [^] for details). The action of this model is 



s = J dY (^$+(y) (A _ A)ci>(y) + 



(3.6) 



Gap {^+{Y) ^\Y) + ^Y) ($+)2(y)) - cl>+(y) Ttr{Y)6{Y - Yo)) 

where G^p is the toy-model vertex that describes the one Pomeron to two Pomerons transition. 
Using the Stratonovich transformation we can linearize the action given by Eq. (|3.6|), namely. 



exp (Ggp Ay ($+)2 ) ~ J df], exp -^r]f - ./2G^^^~AY t],^ 



(3.7) 
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Substituting Eq. ( |3.7|) into Eq. (|3.2| ) and integrating explicitly over $^ we obtain 



(3i 



The (5-function means that the only $ which satisfy to the following (Langevin) equation 

d$(y) = A $(y) - Gsp^'^dY + V2 Ggp (iVF(y) (3.9) 

(where VF(y) is a Wiener process contribute to the functional integral. Eq. (|3.9| ) is actually the 
Langevin equation for directed percolation (see Ref. pH ] for details) which can be re-written in the form 

d $(y) 



dY 



A^Y) - GsP ^'(Y) + C{y) 



(3.10) 



where C is a density ($) dependent Gaussian noise field |2^, which is defined as 

< |C(y)| > = 0; < |C(y) am > = 2G3pHY)6{Y-Y') . (3.11) 



The general functional of Eq. (|3.2| ), we believe, can be reduced to Eq. ( p.lO ). However, we found a 
way how to do this only for Eq. only with addinional sipsimplificationse assume (i) that Eq. ( ^^261) 
has the form 

Li2 G {xi,X2;Y\x'i,x'2;Y') = G {xi,X2;Y\x[,X2]Y') (3.12) 

or, in other words, we consider A{n,u) = A{n = 0, z^) = A(n = 0, = 0) (see Eq. ( |2.1l| )) , which gives 
the leading contribution for the BFKL calculus; and (ii) that 

Xl + X2 



> Xl - X2 ; 



(3.13) 



Indeed, we expect that the typical size of the dipoles will be of the order of 1/Qs{x) {Qs{x) is the satrusat- 
uratione) while the typical impact parameter of the scattering dipole should be much larger (at least of 
the order of the size of the hadron). 

Using these two assumptions we can rewrite the functional integral in the form of Eq. ( 2.15| ) but with 
action Sj that can be written as follows 

d'^xi d^X2 d'^xs 



Si 



2 7ra| 

Nr. 



dY 



111 

/y*^ ry*^ ry*^ 

■^12 -^23 -^13 



{<^>{xi,X2;Y) ■ <^+{xi,xs;Y)^+{x3,X2;Y) + h.c.} (3.14) 



d_ 
dY 



d 



and Eq. ( |2.22| ) reads as 

G^^{xi, X2;Y\x'i, X2;Y') 
Collecting Eq. (3.14), Eq. ( |3.15| ) and Eq. ( 2.23[ ) one can see that the action has a simple form 
So + Si = 



(3.15) 



d"^ Xl d"^ X2 ^+ . d<^ixi,X2;Y) 
-J {xi,X2;Y) — h 



X 



12 



dY 



+ 



dY 



d'^xi d'^X2 d'^xs fas 



2 2 2 
^12 -^23 -^13 



$+(xi,X2;y) {<^ixi,X2;Y) - $(xi,X3;y) - <i> ix3,X2;Y)} 
Ztt 



2 71" Q-2 

+ {^{xi,X2;Y) ■ <^>+{xi,X3;Y)'^+{x3,X2;Y) + h.c.} 



(3.16) 
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Eq. ( 3.1£| ) can be reduce to more elegant form using Eq. ( 3.13| ) and going toomentum representation, 
namely, 

$(xi,X2;y) = x?2 y" (fke'^-^^^^{k,b;Y) (3.17) 

where 6 = (xi + X2)/2 ^ {xi+x^)/2 = {x2 + x^)/2 (see Eq. (^Isj) ). 
In terms of fields ^ (k, b; Y) and <I>^ {k, b; Y) the action looks as follows 

S = So + Si = j cPbcPk^+{b,k-Y) (3.18) 

^d^i^pn _ g I cPk'K{k,k')^{b,k'-Y) + ^.{<I>+6,A:;y)<I>+(6,fc;y) + h.c.]^ 
where K (k, k') is the BFKL kernel in the momentum representation, namely. 

The action S of Eq. ( p.l8 ) looks very similar to the action of the toy model (see Eq. (^^)) and can be 
easily reduced to the Langevin equation for directed percolation 

d^{b,k;Y) _ as f ^2 ,j /^^ u\ ^ !i, 2 7ra? 

dY 

with 



— f ^ — z 

^ / <fk'K{k,k')^{b,k';Y) . ^h,k-Y)^%k-Y) + Q{b,k-Y) (3.20) 



<\ab,k-y)\>= 0; <Ub,k-Y)ah\k'-Y')\>= ^^.^b,k;Y)S^^\b-b')6^^\k-k')5{Y-Y') 

(3.21) 

Actually, we have shown (see Refs.L3,L4) and will discuss below that in QCD we expect a different 
form for the noise , namely, 

< \C{b,k-Y)\ > = 0; (3.22) 
< \C{b,k;Y)C{b',k';Y')\ > = ■ <^{b, k;Y) {1 - <l>{b, k;Y)) 6'^^^ {b - b') 6^^^ {k - k') 6 {Y - Y') 

which belongs to a diffrent universality class than Eq. ( |3.2lD ( ||3l| , |3^ ) ^ . 

The Langevin equation is the one of many ways to describe a diffusion process and the condidereable 
progress has been achieved in this approach, (see Refs.|31, 32, 33, 22, |36|). 

However, we prefer a different way for description of the BFKL Pomeron interactions, which will also 
lead to diffusion equation: the so called generating functional approach. The advantage of the generating 
functional approach is its tansparent relation to the partonic wave function of the fast hadron (dipole) and, 
because of this, this approach can be a source of the new Monte Carlo approach for high energy scattering. 
In this approach we see in the most explicit way our main theoretical problem: the BFKL Pomeron 
calculus provides the amplitude that satisfies the t-channel unitarity while the s-channel unitarity is still a 
problem in the BFKL Pomeron calculus. However, the probabilistic interpretation in the framework of the 
generating functional leads to the correctly normalized partonic wave function which takes into account 
the main properties of the s-channel unitarity as well. 
^Eq. (3.22) for the noise was also suggested in Ref. |3q] 
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3.2 Generating functional: general approach 

In this subsection we discuss tlie main equations of the BFKL Pomeron Calculus in the formalism of the 
generating functional, which we consider as the most appropriate technique for the probabilistic interpre- 
tation of this approach to high energy scattering in QCD. 



To begin with let us write down the definition of the generating functional |17| 

Z{Y - Fo; M) = (3.23) 

/n 
Pn{Y - Yq; xi,yi;...;xi,yi;...;xn,yn) u{xi,yi) <f Xid? yi 

n=L i=l 

where u{xi,yi) = Ui is an arbitrary function of Xi and y^. The coordinates {xi,yi) describe the colorless 
pair of gluons or a dipole. Pn is a probability density to find n dipoles with the size Xi — yi, and with 
impact parameter [xi + yi)/2. Directly from the physical meaning of P„ and definition in Eq. ( 3.23| ) it 
follows that the functional (Eq. ( ^.23 )) obeys the condition 

Z{Y - Yo; [u = l]) = 1. (3.24) 

The physical meaning of this equation is that the sum over all probabilities is equal to unity. 

Introducing vertices for the dipole process: 1^2 (yi_>2(x,y — > xi,yi +X2,y2)), 2^1 (y2^i{xi,yi + 
X2,y2 — > X, y)) and 2 — > 3 (V2^3(xi, yi + X2,y2 — > x'i,y'i + X2, y'2 + 2^3, 2/3) we can write a typical birth-death 
equation in the form 

d Pn^'i Xj, yi] Xji, y-n) /„ 

dY " ^'^■^^> 

= ^ Vi^2^{Pn-i{y;---\Xi,yi;...;Xn,yn)-Pn{Y;...;xi,yi]...]Xn,yn)) (3.26) 

i 

+ ^y2^i {Pn+i{Y;...]Xk,yk;---]Xi,yi-... : x„,y„) - P„(y; . . . ; x^, y^; . . . ; )) (3.27) 

i>k 

+ ^ V2-^3(^{Pn-i{y;xi,yi]...;xk,yk;---;xn,yn) - Pn{Y]Xi,yi;...;xk,yk;---;xn,yn)) (3.28) 



i>k 



Eq. (3.25) is the typical Markov's chain and the fact that we have the correct normalized partonic wave 
function is written in Eq. ( 3.251 ) by introducing for each microscopic (dipole) process two terms (see 
Eq. ( 3.26| ),Eq. ( [3.27] ) and Eq. ( 3.28| )): the emission of dipoles (positive birth term ) and their recombination 
(negative death term). 

Multiplying this equation by the product HILi ^« integrating over all Xi and yi, we obtain the 
following linear equation for the generating functional: 

dZ {Y-Yo;[u]) 



dY 

with 



X[u] Z {Y -Yo;[u]) (3.29) 



XM = j dUd%d% (Vi^2iq^qi+q2) {-u{q) + u{qi)u{q2)) j^- (3.30) 
-V2^l{qi+q2^q) {u{q,) u{q2) - u{q)) ^ ^^^^ ) ; (3-31) 
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Trying to make out presentation more transparent we omitted in Eq. ( |3.3[l| ) the term that corresponds 
to the 2 — > 3 transition (see Ref. for full presentation). We also use notations qi for {xi,yi) and df^qi 
for (Pxi d? yi where Xi and Ui are positions of quark (antiquark) for the colourless dipole. 

Eq. (|3.29| ) is a typical diffusion equation or Fokker-Planck equation with the diffusion coefficient 
which depends on u. This is the master equation of our approach, and the goal of this paper is to find 
the correspondence of this equation with the BFKL Pomeron Calculus and the asymptotic solution to 
this equation. In spite of the fact that this is a functional equation we intuitively feel that this equation 
could be useful since we can develop a direct method for its solution and, on the other hand, there exist 
many studies of such an equation in the literature ( see for example Ref. [|l3|) as well as some physical 
realizations in statistical physics. The intimate relation between the Fokker-Planck equation, and the high 
energy asymptotic was first pointed out by Weigert |^] in JIMWLK approach [^, and has been discussed 
in Refs. ||, p|, |2|. 



The scattering amplitude can be defined as a functional [g^, 19] 

^(X;bi]) = - ^ j ln{xi,yi]..r,Xn,yn]yo)\\-7—Z{Y,[ui])\u,=i(f'Xi(fyi (3.32) 

n=l '' 1=1 ^* 

oo « 

^(-1)" / 7n(xi , yi ; . . . ; Xn ■, y-n ! 



n=l 



The physical meaning of functions 7^ is the imaginary part of the amplitude of interaction of n-dipoles 
with the target at low energies. All these functions should be taken from the non-pertubative QCD input. 
However, in Refs. [l^, ^ it was shown that we can introduce the amplitude of interaction of n-dipoles 
(7n (2^1, 2/1; . . .;Xn,yn',Y) at high energies (large values of rapidity Y) and Eq. ( p.23| ), Eq. ( |3.29| ) and Eq. 
( |3.32| ) can be rewritten as a chain set of equation for 'yn{xi,yi; . . . ; Xn,yn',Y). The equation has the form^ 



djn (qi ■■■,qn) 



dY 

1=1 



n p 

2J2 J d\'d\Vi^2{qu q, q')ln{. ■ -q ■ ■ ) 



ft n It _L n 

^ / d'^q[d'^q2Vi^2{quqi,q2)ln{---,qi---) - ^ j d^qd'^q (3-33) 
i=i i=i 

n „ 

Vi^2{qi; q, q')in+i (. . .g. . .g') - ^ / d^gV2^i (%, qj; q) 

i>j 

n „ 

7„_i {qi...qj ...q) + 2^ d'^qd'^qV2^i {q, qf, q) 7n-i {...qi-.-q) + 

i=i •' 

^ / d'^q V2-.1 {qi,qk;q) -fn{- ■ - qi - ■ - qk- ■ •) 



i,k,i>k 



^This equation is Eq. (2.19) of Ref. but, hopefully, without misprints , part of which has been noticed in Ref. [ p7[ . 
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Comparing this equation for 71 = T^^^ and 72 = T*-^-* one can see that 



V1-.2 



as_ 
2ti 



as_ 
2tt 



K{x,y;z) 



(3.34) 



+ 



2 a| as 
(fx (fy 



T2-,i{xi,yi + X2,y2 x,y) 
{x - yf 

r2_i(a;i,yi + X2,y2 x,y)- 



+ 



(x - yY 

■ (<5(2)(xi-x)<5(2)(yi-y) + 5^^\x2 - x)5^^\y2 - y))) (3.35) 
with T2-*i is given by Eq. (|2.37| ). 

3.3 A toy model: Pomeron interaction and probabilistic interpretation 

In this section we consider the simple toy model in which the probabilities to find n-dipoles do not depend 
on the size of dipoles 18, 20, |2^. In this model the master equation ( |3.29| ) has a simple form 



dZ dZ 

- = -r(i^2)u(i-,.)- + r(2^i)„(i-u)^^^j, 



(3.36) 



and this equation generates: the Pomeron splitting Gp^2P = ~r(l 2); Pomerons merging G2P^p = 
r(2 — > 1) and also the two Pomeron scattering G2P->2P = —T{2 — > 1). It is easy to see that neglecting 



iP' d'^ Z j (duY term in Eq. ( 3.36 ) we cannot provide a correct sign for Pomerons merging G2P^p- 

The description given by Eq. ( |3.36| ) is equivalent to the path integral of Eq. (^^) . To see this we need 
to notice that the general solution of Eq. ( 3.36| ) has a form 

eH(.u)(Y-Yo)z^Y^.u) 



with operator H defined as 

H{u) 

and 



Z{Y;u) 



.r(1^2)n(l-.)|-+r(2^1).(l-.)^ 



Z{Yq-u) 



Introducing operators of creation (a"*") and annihilation (a) 

d 



du 



a = u that satisfy [a, a ] = 1 at fixed Y 
one can see that operator 7i has the form 

n = -r(1^2)a+(l-a+)a + r(2 ^ 1) a+ (1 - a+) 
and the initial state sX,Y = Yq\s defined as 

1^0 > = e^"-(^''"i)|0 > 



(3.37) 

(3.38) 
(3.39) 

(3.40) 
(3.41) 
(3.42) 



1 o 



with the vacuum defined as a|0 >= 0. 

We need to discretize the development operator of Eq. ( 3.37 ) with TL given by Eq. ( 3.41| ), namely, 



N 

^n{Y-Y,) = e^^^ . . . e^^^ = n (1 + W^^) (3-43) 

and introduce the coherent states [^] for a certain intermediate rapidity Y^- = Yq + j Al" as 

> = S'^i'-^-'i', |o > (3.44) 



where are arbitrary complex numbers. The initial state of Eq. ( |3.42 ) can be written as 



\4>^{J^)>=\Ttr> ■ (3.45) 
The unit operator in terms of the coherent states can be expressed as 

1 = f^e-^^<^l + <^^+<^.*)|<A,><<A,| (3.46) 

We want to calculate matrix element of some operator A between states of initial Yq and final Y 
rapidity < y|A|yo >• This can be written as 

< y|^|yo >=< ( / ^^^e-<^^^^+<^^+'^?-) |0y X </.y|| 1^0 > (3.47) 

here we denote |y >= \4>y >. Next we use the development operator given in Eq. ( |3.43| ) to find < Y\Yq >. 
We split the rapidity Y — Yq to N intervals. And insert the development Eq. ( p. 431 ) and unit Eq. ( |3.46D 
operator between the states of intermediate rapidity 

N 

<Y\Yl {1 + HAY) \Yo> (3.48) 

We look at 

< {1 + n AY) > = exp + + ^j+i - - + 

X (1 + n{<P*+„(P,)AY) 

= exp {-0*+i(0j+i - 0j) + - (Pj} (1 + ?^ (</.*+!, 0^) AY) 

= exp {-0*+i(0j+i - 0j) + - exp (W(0*+i, cPj) AY) (3.49) 

Now we redefine an arbitrary function (pj as 

= $+ = 0*-l (3.50) 
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and rewrite Eq. ( |3.49| ) in terms of and <I>^ 



exp {l>++i($,+i - ^>,) + + 1, -1>,)} 

''^$++i('&i+i-'J>,) 



(3.51) 



exp 



Ay 



N 



Summing over all rapidity intervals we have 

< y|^|yo > ~ n / d^^d^jA{Y)e^ 
where A(Y) is the expectation value of the operator A at the final rapidity Y, and 



(3.52) 



(3.53) 



In the continuous limit this becomes 



< Y\A\Yo > 



J V<^+V<^> A{Y) e" 



with 



— r(i ^ 2)$+$ + r(i ^ 2)$+($)2 + r(2 ^ i){^+f^ - r(2 ^ i){^+f {^f ) dY 



(3.54) 



(3.55) 



This action is almost the action of Eq. ( p.6| ) for r(l — > 2) = r(2 ^1). In the toy-model the difference 
between these two vertices is the normalization problem of function and ^. In our approach they 
are normalized in the way which allows us to treat them as probabilities (see Eq. ( |2.19| ) and Eq. (2.20)). 
However, Eq. (3.55) includes the new interaction: the transition of two Pomerons to two Pomerons. The 
sign is such that this interaction provides the stability of the potential energy. Indeed this term is responsible 
for the increase of the potential energy at large values of both and ^. 

Comparing Eq. ( 3.55| ) with Eq. ( p.l| ) one can see that we build the partition function and the thermo- 
dynamic potential using the generating functional. It means that our Eq. ( p. 29 ) is equivalent to statistical 
description of the system of dipoles. 

Eq. ( 3.36| ) is the diffusion with the u dependence in diffusion coefficient. In terms of the Langevin 
equation Eq. ( |3.55 ) generates the noise term of the Eq. ( 3.22| ) type. 

To out taste this equation is simpler than the Langevin equation of Eq. ( |3.9| ) and it will be easily 
generalized for the case of QCD. For u < 1 the diffusion coefficient is positive and the equation has a 
reasonable solution. If n > 1, the sign of this coefficient changes and the equation gives a solution which 
increases with Y and Z(Y) cannot be treated as the generating function for the probabilities to find n 



an 



dipoles (Pomerons) (see Refs. [10, for details). The same features we can see in the asymptotic 

solution that is the solution to Eq. ( 3.36| ) with the l.h.c. equal to zero. It is easy to see that this solution 
has the form 



Z{u;Y^ oo) = ^— ^ ; with . = = ^ » 1 (3-56) 

One can see that for negative k this solution leads to Z > 1 for u < 1. This shows that we cannot give a 
probabilistic interpretation for such a solution. 



4. Probabilistic interpretation in QCD 
4.1 Several general remarks 



Eq. ( 3.34 ) has a very simple physical meaning describing the Pomeron splitting as the decay process of 
one dipole into two dipoles. It turns out that the vertices for 2P 3P and 2P 2P can be easily 
understood as a dipole 'swing'. What we mean is that with some probability two quarks of a pair of 
dipoles can exchange their antiquarks to form another pair of dipoles ||2^. Naturally, this process has 
suppression and it correctly reproduces the splitting and rescattering of two Pomeron that has been 
explicitly calculated from the diagrams ||^]. 

Eq. {^MY 

is more difficult to view as the vertex for the transition of two dipoles into one. Indeed, the 
integral over coordinates of the produced dipole is positive, namely / V{2 —>■ l)cfx(Py = f T{2 —>■ l)(Px(Py. 
However, Eq. ( p.35| ) generally speaking leads to a negative vertex in some regions of the phase space. Here, 
we want to point out that the key problem is not in the probabilistic interpretation of the microscopic 
process of two dipole to one dipole transition but the fact that a negative vertex V{2 — > 1) means that 
in some kinematic region we have a negative diffusion coefficient which results in a solution that increases 
at large values of rapidity Y. To save such theory we need to introduce other Pomeron interactions like 
2P 3P and/or 2P 2P transitions. 

In Ref. H, H] the attempts were made to deal with such Pomeron interactions. It turns out that at 
large values of Y the 2 — > 1 process contribute in a very limited part of the kinematic region with very 
specific function u{x,y). In this particular region the vertex V{2 — > 1) is positive. Therefore, we could 
use the probabilistic interpretation but we need to study this process better and deeper to obtain the final 
result. 

In the toy model it has been shown | p!o| , [l^ that we can generate the 2P — > 3P vertex without the 
process of two dipoles to one dipole transition. The transition of two dipoles to two or more dipoles also 



leads to this vertex. The similar ideas are developed in Ref. |45| in QCD. However, we need to pay a price: 



the contribution to Eq. ( ^.29 ) will be negative to the correct sign for 2P P interaction. In other word, 



we can add to Eq. ( ^.291) the contribution 
» 2 2 

/ (fxi (fyi Y{ d^x'i Sy\ V2-.i{xx,y\ + Xi,y2 x'i,y[ + x's, ^2) 

i=l i=l 

^This equation is quite different from the equation which is obtained in Ref. (see also Ref. p7[]) . The main difference 
stems from the correct use of the BFKL Pomeron Calculus for determining this vertex while in Refs. ||2C|, [27| the Born diagram 
was used which does not and cannot give a correct expression. 




Figure 4: The diagrams for 2P to 2P process. Fig. ^a shows the diagram which is of the order as/N^ while Fig. ^a 
leads to a contribution oc a^. Fig. ^c shows the first diagram for the 2 dipole to 2 dipole amplitude which square 
gives the probability for 2P to 2P process. 

Id d 

u{x'^,y[)u{x'^,y'^)-— -— (4.1) 

2du{xi,yi) du{x2,y2) 

which will give the 2P P vertex in the form 

r{2P^ P){xi,yi + X2,y2^x[,y[) = (4.2) 

- j d^x2 (fy'2 ^2^2(2:1, yi + X2,y2^ x[,y[ + x'a, ya) 

Therefore, we have to assume that V2->2 < 0. As we have discussed, generally speaking, it means that 
our problem has no solution. However, in QCD the situation is much better. Indeed, to generate a correct 
2P — > P vertex we need to introduce V2-+2 oc a'g. The Feynman diagrams for this transition is obvious 
(see Fig. (Q-b). However, the main contribution stems from the diagrams of Fig. (^a) - type which are of 
the order of as/N^- Therefore, diagrams of Fig. ^b - type are small corrections to the main contribution 
and could be negative, in spite of the fact that the diagrams shown in Fig. (^b) actually gives a positive 
contribution. 

The situation with ¥2^1 is actually more dramatic since even in this paper we have discussed four 
different expression for this vertex: (i) Eq. ( ^.35 ); (ii) Eq. ( |2.36| ) leads to the vertex of the following form 



(see Refs. |22, 23, 20] where this form of the vertex has been discussed in details) 

T2-,i{xi,yi;x2,y2\x,y) =V^V^ j (fzK{x,y;z)Go{xi,yi\x,z)Go{xi,yi\z,y); (4.3) 
(iii) Eq. (|1|) leads to 

T2^i{x,z;z,y\x,y) = K {x,y\z) (4.4) 

and the expression for vertex V2^i one can obtain substituting Eq. (|4|) into Eq. ( ^35)) ; (iv) Eq. ( p^ 
suggests the simple vertex 

V2^i{b,k;b,k\b,k) = (^^) (4.5) 
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These expressions are equivalent in the framework of the BFKL Pomeron calculus but lead to quite different 

|) does not generate 2 



physics. For example, Eq. ( |4.3| ) does not generate 2 — > 2 dipole term in the functional integral [27| and, 
therefore, the noise term has a form of Eq. ( p.21 ), while all other expression correspond to the noise term 



of Eq. (3.22) which belongs to a different universality class. 

As we have mentioned that in the kinematic region where L{Y) ~ 1 (see Eq. (p. 7])) the 2P 2P 
transition] gives only neglitransitionibution. However, it changes crucially the behaviour of the scattering 
amplitude at ultra high energy. 

Our strategy of the discussion looks as follows. First, we are going to discuss the two Pomeron to 
one Pomeron transition directly in terms of the dipole approach calculating Feyman diagrams. The key 
question that we want to answer is what process of the dipole recombination is described by the diagram 
that contribute to the 2P P vertex. Second, we wish to find what term in x['u] (see Eq. ( 3.30D ) is 
responsible for this contribution. Third, we find the solution for the toy model that includes this term. 
Finally, we develop the generating functional approach for the approach that correspond to the action 
given by Eq. ( pl|) which we think will be a practical way in searching the solution. 



4.2 Scattering of two dipoles in QCD 

We start to approach the problem of taking into account the two Pomerons to two or more Pomerons 
transition with the clear understanding of 2P — > P merging in the dipole approach. In this section we will 
obtain the same Eq. ( p. 37 ) for r2^i but directly from the dipole picture of interaction without using the 



trick suggested in Ref. |22|. Namely, we will calculate the contribution of the simple diagrams like that of 
Fig. 18-a but directly from the dipole picture of interaction [17| shown in Fig. [l^b. 



Starting to calculate the diagram of Fig. 18 we would like to draw your attention to the fact that 
these diagrams actually describe three different processes. In the Born approximation all of them look 
quite similar but if you imagine, that gluon with momentum k and I can produce more gluons (see dotted 
lines in Fig. Fig. 19 and Fig. 2C) , we can see a great difference between three diagrams. Indeed, the 



diagram of Fig. 18 represents the process where two bunches of additional gluons could be produced from 
gluons with momenta k± and l±. This diagrams give a positive contribution and has clear probabilistic 
interpretation. Fig. ^ shows the process where none of additional gluons can be produced. The gluons, 
that were produced in the process of the diagrams of Fig. |l8|, contribute to the structure of the BFKL 



Pomerons. The diagram of Fig. 19 corresponds to a process where only one of the bunches of gluons is 
produced while the second one is absorbed by the BFKL Pomeron. This diagram belongs to an interference 
type of the diagram which can have a negative sign. 

It is clear that these extra diagrams correspond to different AGK cuts of the BFKL Pomeron diagrams. 
Fortunately, they differs by the factor and sign in front of the same expression [BTl, namely. 



Fig.^ : Fig.\l^ : Fig.\]^ = 1 



(4.6) 



Before approaching our problem we want to show how one can calculate well known Lipatov vertex 
using light cone perturbation theory. The calculations are performed in the infinite momentum frame, 
large Nc, eikonal and Regge limits (E^, E3|. 
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Figure 5: The diagram for the gluon exchange in the infinite momentum frame. The first diagram survives while 
the second one gives a smaU, neghgible contribution at high energies. 

We re-calculate the simplest diagram of one gluon exchange in the infinite momentum frame technique 
to illustrate the simple space-time structure of the high energy scattering that first have been formulated 



by Feynman and Gribov |42, |44[ . 



It is well know that the one gluon exchange in Feynman diagram approach is equal to 

Ag = ^ (4.7) 

In the infinite momentum frame we need to sum the contributions of two diagrams shown in Fig. ^ which 
are equal to 



2 Eq Ep^ + Ep^ — Ep2 — Epi^ — Eg 2Eq _ pi,± _ ^ 1± ' 

2Eq Ep^ + Ep^ — Epi_^ — Ep^ — Eq 2Eq 2Eq ( / ) 1 ( ) 

The difference between these two diagrams is the following. The first one describes the process with 
natural ordering in time. Indeed, the typical time for the relativistic particles is E/p'^ where E is the energy 
of the particle and p± is its transverse momentum. Therefore, the first diagram describes the process of 
the emission of long living gluon with the lifetime Tq oc Eq/q'j_ by the parent particle with its lifetime 

oc Ep-^/p1j_ ^ Tq. In the second diagram, the particle with short lifetime oc Ep^/p'^_^ <C emits 
the gluon with larger typical time. For emission of large number of gluons we have a space-time picture 
shown in Fig. |6|. The square of these ladder-time diagram we will call the BFKL Pomeron. The general 
feature of this diagrams that the cross section can be written as the product of two factors: the square of 
the multi-gluon wave function {^) and the cross section of the interaction of the slowest, 'wee' gluon with 
the target 

a = \^^^a{Gluon- target) (4.10) 
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time 



Figure 6: The space-time structure of tlie BFKL Pomeron in QCD. 




Figure 7: The total cross section of two quarks on one scattering 

where stands for integration over all kinematic variables of the wee gluon and is the wavefunction 
integrated over all kinematic variables of the gluons that are faster than the wee one. 

Before starting our calculations we want to show that this kind of diagrams indeed give the leading 
logarithm contribution. We consider the total cross section shown in Fig. 0. 

We want to perform the integration over plus and minus components of all internal momenta leaving 
the transverse integration. It is convenient to do the integration in terms of Sudakov parameters on the 
complex plain. We assume the Regge limit, PiPs^ ^ ^2^3^ = s/2 and forward scattering, i.e. Q± = 0. 
Thus the internal momenta can be written as 



nr 



Pi ;Prm a,c ; p.-m-q 




Figure 10: A diagram contributing to the effective Figure 11: A diagram contributing to the effective 
vertex. vertex. 



m'' = amp't + f3n,P^ + ml (4.11) 
The total cross section have denominators that come from gluon propagators 

= akPkS - kl (4.12) 

= ttmPmS - m]_ (4.13) 

l^ = ai[5is-ll (4.14) 

{k - mf = {uk - am){l3k - Pm)s -{k- m)l (4.15) 

(^l - Qf = (ai - agm - Pq)s - ll (4.16) 

(k + Qf = {ak + aQ){Pk + f3Q)s - kl (4.17) 

{k-m + Qf = {ak -am + aQ){Pk - Pm + Pq)s - {k - m)i (4.18) 



and the 5-functions (from the cuts) of 



{P2 - If = -Pis - ll (4.19) 
{pi - kf = -(3ks - kl (4.20) 
{P3 + mf = OmS - m]_ (4.21) 

In forward scattering the transverse part of the transferred momenta is zero {Q±_ = ), from the on- 
shellness of the outgoing quarks it can be seen that aq and /3q can be neglected in further calculations. 
According to the Cutkosky rule the amphtude should be squared and integrated over ah internal momenta. 
If the integration is performed in terms of Sudakov parameters it reads 

(fk dH (fm = {^'^dak djJk dai dPi dam d(3m d^k± d^l± d^m± (4.22) 

The integration over transverse momenta do not bring the leading logarithm of energy and thus may be 
omitted in our discussion of the origin of the In s contribution. 

First we do the integration with respect to Pk which gives Pk = ~kj_/s and factor —l/s from 5^((pi — 
kf) (see Eq. 4.19| ). Next the integration over Pi is performed, this gives at Pi = —l\/s and factor —l/s 



from 5'^{{p2 — if)- The (^-function of (ps + mf gives am = +m^/s and factor l/s. For t-channel gluons 
the propagators can be written as follows 



e = akPkS -kl= ak{^)s - fei ~ -kl (4.23) 

w? = amPmS -mj_ = —^PmS - mi ~ -m\ (4.24) 

f = aiPis -ll = «/(—)« -ll^ -ll (4.25) 

s 

Similarly one finds that {k + Qf ~ -kl, (/ - Qf ~ -ll . 

After this we proceed to calculation of other denominators are given by 



{k - mf = (ak - am){Pk - Pm)s - {k - m)i = (4.26) 



ak )( — - - Pm)s - {k- m)_L = akPm - {k - my 

s s 



and the integration with respect to ak is done on the complex plain adding to denominator of Eq. 4.26 
some small ie. The contour integration gives a residue at ak = {{k — m)l —i(\jPmS and the factor l/[/3ms]^ 
before the amplitude (/9^ squared is because of the propagator of gluon m — k — Q). Next, {k — m + /)^ 
gives ai = [{k — m + 1)1 — ie]/Pm and another factor if 1/PmS. 

We are left with some function of only transverse variables times dPm/Pm- '^^ notice that the two 



triple gluon vertices in the center of Fig. ^ result into Pm in the numerator (see Eq. 4.44| ) canceling two 
powers of Pm, than we end up with dPm/Pm which brings the logarithm of the energy. 



We start with three types of diagrams of quark-quark scattering contributing to the effective vertex 



(see Fig. g- 11). This result will be useful for our further considerations. The amplitude depicted of Fig. 
S reads as 



^3 ujPi) ^M^b ^(Pi - -g^u 1 u{pi - m) _ u{pi -m-q) u{p^) ^^^^ 1 u{pz + m) 



where Di and D2 are given by 
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Pi- + P3- 



{Pl 



m) 



PU 
{Pi 



m 
^2 



771 ]_ 



P3- 



{P3 + ITT-Y 



PU 



rriA 



P3+ + "T-H 



(4.28) 



{pi-m- q)]_ (p3 - m)^ (pi - m)\ {ps + m)^ 



+ 



+ 



Pi+ -m+-q+ P3+ - m+ pi+ - m+ ^3+ + m+ g+ 
In estimates of D2 we used the fact that initial light cone energy equals final light cone energy at high 



energies. Summing Eq.(4.27) over quark polarizations (see Ref. |^3|) we have contribution only from plus 
and minus components of 7^^ and 7^ respectively. Namely, 



u{p) n u{p-m) 



(4.29) 



where is defined such that n^kn = for an arbitrary 4- vector k^. Now we can rewrite Eq.4.27 as 



-m+ —1 
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ml 



Here we used the transverse condition 

1 . 







(4.30) 



(4.31) 



One should consider another diagram depicted in Fig. but it can be readily seen that both of its 
denominators are of the order of Di ~ D2 ~ m_ = / . Thus it is has an additional m_ suppression 
comparing to the amplitude of Fig. ^ which will not give the logarithmic contribution oc Y = Ins. Next 
we consider a diagram shown in Fig. Il^ given by 



9 



uipi) a^c uiPl -q) 1 u{pi 



/Pi+ VPi+ - 1+ Di ^/pi- 
where Di and D2 are given by 



.M^b ujpi - m-q) u{p2) ^^^^f, u{p2 +m) 1 



Vpu 



/P3- 



^P3_ + m_ D2 



(4.32) 



Di = pi_ + P2- 

D2 = Pi- + P2- 



m ■ 
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{Pi 



Q+ 

\2 
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(4.33) 
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Summing over quark polarizations we get the resulting expression for the amplitude shown in Fig. 10 

t^t''0tVn+n_ ~f 6- = -tH'^<g>t'g'2^-^e. = -th'^t'gH^^ (4.34) 

qj^ m]_ m+ m]_ m]_ qj^ 

We see that the contribution of Fig. IC is equal to that of Fig. |8| but opposite in sign. We also notice that 
they enter with different color matrices f^t^ and t^t"' (8) respectively. Using the basic relation for Lie 
group generators = if^^^t^ we can write the contribution from the two diagrams as 



The last diagram to be considered is shown in Fig. |Tl and given by 



where Y'^'^^''^ is a regular triple gluon QCD vertex. 



7' 



(4.35) 



(4.36) 




The light cone energy denominators Di and D2 are readily 
calculated 



Di = pi_ +P2- 
D2 = Pl- + P2- 



fa-fc)i ^Jl 

(Pi-fc)i gi _mi ^_!!i(4.37) 

pi+ - k+ q+ m+ m+ 



Figure 12: Schematic representation of 
the effective vertex. 



The triple gluon vertex can be reduced as follows (see IQ) 

= gig^'^^'ik + mf + g^'''{-m + qY' + /"(-g - A;)'''] 
= [/^'2r + g^'^^-m + qY' - /'"2F'] ~ /^'2/t'"(4.38) 



With these simplifications Eq. 4.36 can be written as follows 



-k+-m+ -1 -1^,, 



k\ m\ m+ /c+ 



1 fc_Le_L 



ml 



(4.39) 



The resulting contribution from diagrams Figs. p[-p!l| is schematically shown in Fig. 12 and reproduces 
the effective Lipatov vertex 

(4.40) 



abc^c ^ ^6„3^ J_ ( k±e± _ q±e±\ 
2 ^ A;2 



J 



We continue with the light cone perturbation theory and consider a second quark emitting gluon I 
that hits gluon q, as shown in Fig. |l^. It is worth mentioning that the gluons / and q should be coupled 
to the quark 1 and quark 2 lines via the effective vertex found before. For simplicity, we consider only one 
type of the diagrams contributing in this case and generalize the result to all the possible couplings. There 
exists a lot of ways to perform the light cone energy cut when one adds second scattering quark. As it was 



shown before that cuts double the same lower gluon are suppressed in Regge limit. The amplitude shown 
in Fig. 13 is written as follows 



U{pi) u u{pi-k) -g^^t 1 u{p2) „ U{p2-l) -9(3(3' 1 

T — — — — 

^/pT+ y/pi+ - k+ k+ Di ^/p^ Y/p2+ - 1+ h D2 



D3 ^/P3+ m+ 
where the denominators D1-D4 are given by 



T.^-M-a 1 ujPs) ^ -g^,^^' u{p3-m) -g„„, 1 /^/^ 

~ ' / n ^ 
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(4.41) 
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(4.42) 



In Eq. 4.41 we omitted coupling constants and color factors that will be easily restored at the end of the 




Figure 13: One of the diagrams contributing to 2P ^ P vertex calculated in light cone perturbation theory 



calculations. The amplitude Eq. 4.41 can be further simplified as follows 

-W -k+ -913(3' -h p^Va' -m+ -g^^> -g° 



n 



I, 



q+q- 



(4.43) 



where vector is defined as n^/c^ = for any arbitrary vector k^ . The triple gluon QCD vertex 
is given in Eq. 4.38| . The vertex E'^'^'" in our case can be simplified as follows 



E 



g'^(''{l-qY +g--'{2q + ir +g' 



-q - 2Z)" ~ 5"'" {2q + I) 



Using this we rewrite Eq. 4.43| as 



+ 71^2 



1 1 1 (2g + /)_ , , 
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(4.44) 



(4.45) 



We note a plus component of q in the denominator, this factor is not related to integration in transverse 
coordinates and leads to a logarithm of energy as was shown before for total cross section, thus it can be 
omitted at this stage and safely restored at the end of the calculations. Performing similar transformations 
for an amplitude of two gluons being emitted from quark (antiquark) lines (see Eig. ^ p!o| ) we get 

1 1 1 



+ 2-2 



q\ m" 



-2q± ■ 



with the same color factor as of Eq. 4. 45. Thus, Eig. |1^ gives 
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Ij mi 
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(4.46) 



(4.47) 



One should add another one, where momenta q and I are interchanged, to this expression, then it gives 
an effective vertex for the scattering of two with one quarks (with color factor, powers of g and plus- 
minus components restored). Now we want to find the cross section that corresponds to two dipoles being 
scattered with a target. We use a mixed representation where we Eourier transform all transverse momenta 
to coordinates of dipoles (quarks). We assign transverse coordinates Xi [yi) to each quark (antiquark) line, 
and coordinate z for the emitted gluon. The Eourier transform is performed separately for (fik_i_ and dPq± 
integrations . We start with Eq. |4.45| 

cP'kj^ (Plj_ cPmj_ 



I 



(27r)2 (27r)2 (27r)2 



-ik±yi 



-il I X 



i_X2 



-ilA_y2 



To simplify Eq. 4.4£ we introduce the gluon propagator in the coordinate space 



Dg{x,x) 



d^^l. 1 Jkdx-x') _ ln(|x-xV) 



(4.48) 



(4.49) 



{2ttY kY 27r 
where ji is an arbitrary small mass that we need to send to at the end of our calculations, and the integral 
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2-K {x - x'f 



(4.50) 



which win be interpreted later. In terms of Eq. 4.49 and Eq. 4.5C| we can rewrite Eq. 4.48 as follows 

'{z-xi)°' (z-yi) 
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{z-xiY {z-yiY 
In a similar way we rewrite Eq. [4.46 in the form of 



{Dg{z, x-i) - Dg{z, ys)) {Dg{z, X2) - Dg{z, ^2)) 



(4.51) 
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{Dg{x3,xi) - DG{y3,xi) - Dg{xz, yi) + -Dg(2/3, yi)) X2) - Dg{z, ^2)) 



(4.52) 



Eqs. 4.51- 4.5^ should be multiplied by the wavefunctions of the original dipoles ^^^\xi,yi), their factors 
have transparent physical meaning, namely 



{z-xi)" (z-yiY 



(z-a;i)2 (2;-yi)2_ 

{Dg{z,X3) - DG{z,y3)) 



describes the emission of gluon z by dipole {xi,yi) 



describes the interaction of gluon z with target dipole {x3, ys) 



(4.53) 



It should be noted that Eq. 4.51 should be multiplied by 2 before adding to Eq. 4.52. This could be 
explained as follows, let us go back from large Nc limit to A'^c = 3 and add a color singlet corresponding 
to large Nc limit. We want to consider a combinatorial factor for all possible emissions of two gluons from 
dipole {xi,yi) (see Fig. 14- Fig. ^) fixing the colour in the final states. It is clearly seen from the picture 
that the upper figure that corresponds to Eq. 4.45 , has as much as twice options for two gluons to be 
emitted. 
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Figure 14: Colour structure of the reaction 
q{blue) + q{antiblue) ^ qq + G(q + m) qq + 
G{q) + G(m). B, G, R{B, G, R) stand for blue, green 
and red quarks (antiquarks), respectively. 
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Figure 15: Colour structure of the reaction q{red) + 
qiantired) qq + G{q + m) qq + G{q) + 

G{m).B, G, R{B, G, R) stand for blue, green and red 
quarks (antiquarks), respectively.labelCol2 



Using this fact we may add Eq. 4.51 to Eq. 4.52 
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{DG{z,xi;x3,y3) + L>g(z, yi; X3, ya)) {Dg{z,X2) - DG{z,y2)) (4.54) 
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Figure 16: Colour structure of the reaction 
q{blue) + q{antiblue) — > qq + G{q) + G{m). 
B,G,R{B,GB,G,R{B,G,R) stand for blue, green 
and red quarks (antiquarks), respectively. 



Figure 17: Colour structure of the reaction 
q{green) + qiantigreen) — > qq + G{q) + G{m). 
B,G,R{B,G,R) stand for blue, green and red 
quarks (antiquarks), respectively. 



where Dg{z, xi;xs, 2/3) stands for interaction of dipole (z, xi) with the target. One should also add a term 
similar to Eq. 4.54| with dipole 1 interchanged with dipole 2 {xi,yi <-> X2,y2)- Now we are in position to 
square the summed diagrams. It can be easily shown that all types of interference terms vanish after the 
integration over z. Thus we are left with each term squared and the final expression reads as 



dY' [ (fz{2j''^{x2,y2\z)K{xi,yi\z)[j''^{xuz\x3,y3)+l'''^{yi,z\xs,y3)] + (1 ^ 2)} 



(4.55) 



{x,y\x ,y 



^ 47r^ 



'Go(a;, y') with Go being the initial Green function given by Eq. 2.24 



where 

and K{x,y\z) is given by Eq. p. 34 . J^^{x2,y2\z) describes the interaction of a gluon with coordinate z 

(Xl - 2)2" 



with dipole {x2,y2) (see Fig. |2^). It is equal to 

7^^(xi,yi|z) = 



8 V27r/ 



(yi - 



(4.56) 



Eq. ( 4.55| ) should be multiplied by the factor Ins which is originated from the integration over the energy 
of produced gluon with the coordinate z. 



One can see from Eq. (4.54) that we obtain the kernel of the BFKL equation for the interaction of 



one dipole {xi,yi) in Fig. 18-d with the target dipole {x,y) (see Fig. ttq-a). Indeed, neglecting factor 



Dg{x2, z) — Dg{z, y2) in Eq. ( 4.54 ) we have the following answer taking derivative with respect to F = In s 
from \A{Eq. (^))p 



d\A{Eq. (g:54|) )P 



dY 



(4.57) 



as 
27r 



d'^z\^o{x2,y2)\'^ K {xi,yi\z) (7'^^(3;i, z|x3, ya) + 7^^(^, ^1^3, ya) ) 
In Eq. ( [4.57 ) one can recognize the BFKL equation with the kernel ^ K {xi^yi\z). 



00 



4.3 Contribution to the generating functional 



Eq. ( 4.55| ) gives the expression for the transition (xi,?/i) + {x2,y2) {xi,y2) + {xi,z) + {z,yi) given 
by the diagrams of Fig. 18-c and Fig. 18-d. This diagram gives a positive contribution but the sum of ah 
diagrams (see Fig. |l8|, Fig. 19 and Fig. leads to a negative contribution due to the AGK cutting rules 
(see Eq. (|4.6| )). The advantages of this direct calculations are clear: first, we obtain the vertex for 2 — > 1 



transition in more convenient form than it was calculated before (see Refs. [P2| , 23, |2^). From Eq. ( 4.55 ) 
we see that this vertex is equal to 



r2^i {{xi,yi) + (x2,y2) {x,y)) 
(^sNc ro RA 



27r2 



K {xi,yi;x2,y2\z) 



(4.58) 



2^2 {'^'^ {xi,yi\z) X K {x2,y2;z) + xi ^ X2,yi ^ y2} 
It should be stressed that Eq. ( [4. 581 ) leads to a vertex which is positive in the full phase space. 

The second advantage is that we can understand what process in terms of dipole these diagrams 
describe. Indeed, at initial rapidity Y we have two colour dipoles:(xi, yi) and {x2,y2)- The diagram of 
Fig. 21 -a describes the transition of these two intitial dipoles which goes in two initial 

(xi,yi) + {x2,y2) {xi,yi) + {x2,z) + {z,y2) (xi,y2) + {x2,z) + (z,yi); (4.59) 

and the target dipole {x,y) interacts with dipoles: {x2,z) and {z,y2) during the first stage of the process 
or with {x2,z) and {z,y2) in the last stage. 

Therefore, in terms of the functional of Eq. ( p. 32 ) these diagrams can be written as 



r2_i(^^. ( p8| )) (7(xi,z) + 7iz,yi) + jix2,z) + 7(2,7/2)) 



^l{xi,yi) S-f{x2,y2) 



(4.60) 



The form of Eq. ( [4.6[)| ) suggests that this term can be originated by the following term in the Hamiltonian 
X[u\ (see Eq. and Eq. (|]3|)) : 



1 s s 

-T2^i{Eq. (|4.58|) ) [u{xi, z) u{z,yi) + u{x2, z) u{z,y2)] —, -— r 

2 du{xi,yi) du{x2,y2) 



(4.61) 



The Feynman diagrams that correspond to this term are shown in Fig. ^b and Fig. ^c. However, Eq. (4.61) 
cannot be correct since it leads to a positive two dipoles to two dipole amplitude. Indeed, for two dipoles 
to two dipole amplitude we have the same three diagrams of Fig. 18, Fig. 19 and Fig. 2C types. The sum 
of these diagrams gives the negative sign (see Ref. [39| where this problem was studied in details ). 



Repeating the same calculation which we made, summing Fig. |18|,Fig. |19| and Fig. |20|, one obtain for 
2^2 dipole amplitude the following answer 



1 



T2^iiEq. ( p8|) ) i^(xi,z)h(z,yi) + j{x2,z) + 7(^,^2)) + 



(4.62) 



'y{x2,z){jiz,yi) + 7(2,2/2)) + 7(^,^1)7(^,2/2)} 



^l{xi,yi) Sj{x2,y2) 



O A 






Figure 18: The diagram of the 2 Pomeron into one Pomeron transition and its dipole interpretation: diagonal 
contribution. 



Eq. ( [4.58| ) and Eq. ( [4.62| ) suggest that all dipoles that have been produced in both stages of the 
process (see Eq. ( 4.59D ) interact with the target. It means that the contribution to x[n] (see Eq. ( |3.29| ) and 
Eq. ( |3.3[1| )) looks as follows 



2 



{Eq. ( |4.58| )) u(x2,z)u{z,y2)uixi,y2)u{z,yi] 



6u{xi,yi) 5u{x2,y2) 



(4.63) 



Eq. ( [4.63| ) describes the process of decay of two dipoles into four dipoles^. As one can see from Eq. ( p. 261 ) 
- Eq. (|]2|) each process for the dipole transition generates two terms in Markov's chain which lead in 



Eq. ( |3.26| ) - Eq. ( 3.28 ) to the positive contribution which represents the increase in the number of dipoles 
due to the elementary process ( birth term)and a negative contribution which describes a decrease of the 
number of dipoles as a result of the elementary process (death term). The death term which corresponds 



As far as we know A. Kovner was the first who discussed the possibiUty that two Pomeron to one Pomeron merging is 
related to the decay process in the dipole approach. 



or 
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Figure 19: The diagram of the 2 Pomeron into one Pomeron transition and its dipole interpretation: non-diagonal 
contribution with elastic rescattering of the dipole xi,yi). 



to the birth term given by Eq. ( 4.63| ) has the form 



1 S S 

- / (fz T2^i (Eq. i ^M )) u{xi,yi) u{x2,y2) j—, r t-t r 

2 J 6u{xi,yi) du[x2,y2> 



where 



J <fz T2^i {Eq. (14581) ) 



(4.64) 
(4.65) 



asNc 
27r2 



2(7 (xi,yi|x2) , + 7 (xi,yi|y2)) X / d z K {x2,y2; z) + xi^X2,yi^y2 



It is easy to see that for the case when two dipoles are the same, Eq. ( 4.64 ) describes the diagram of Fig 
Finally, the equation for x[u] (see Eq. ( 3.30| )) has the form 

6 



X[u] = { I d^xd^yd^zK {x,y\z) {-u{x,y) + u{x , z)u{z , y)) 



2tt 



Su{x,y) 



(4.66) 



1=1 



X iu{xi,z)u{z,yi)u{x2,z)u{z,y2) - u{xi,yi) u{x2,y2)) 



6u{xi,yi) Su{x2,y2) 



as 1 



^ 1=1 



d^XiSyid^z {K{xi,y2\z) + K{x2,yi\z)) 




Figure 20: The diagram of the 2 Pomcron mto one Pomeron transition and its dipole interpretation: non-diagonal 
contribution with the interference between the decays of two dipoles. 




Figure 21: The diagrams for the amphtude of two dipole with one dipole interaction with emission of one extra 
gluon. 



X (1 - u{xi, 2/2)) (n(x2, 2/1) - u{x2, z) u{z, yi)) 



6u{xi,yi) 5u{x2,y2) 



where the last term is the transition of two dipoles to three dipoles taken from Ref. [20|. We introduce 
this term for completeness of presentation as well as to demonstrate that a part of terms for the Pomeron 
interactions induced by the second term in Eq. ( ^.66| ), have a contribution which is suppressed by the 
extra power of the QCD coupling in comparison with the same interactions generated by the last term in 
Eq. (3). 



4.4 Generating functional in the toy model 

In this subsection we return to the simple toy model to investigate the specific properties of the generating 



Figure 22: The diagrams for the amplitude which describes the death term in the equation for the generating 
functional in the case when dipoles {xi,yi) and (x2,j/2) are equal. 



functional given by Eq. ( [4.66| ). The equation for the generating functional in this model has the form 

(4.67) 



|_ = _r{1^2)u{l-u)^ + {r(2^1)n2 (1-^x2) +r(2^3)n(l-n)2} 



du 



(du) 



The asymptotic solution to this equation can be found solving the equation with = 0. In semi- 



classical approach for which Z{Y = oo,u) = e**^") with ^u,u *S ^(u) is equal to 



where 



r(i 



Kl 



r(2 



Prom Eq. (|^) 



2) 1 

« — > 1; 
1) ai. 



Z{oo, u) 



r(2 ^ 3) 



K2 



r(2 ^ 1) 



oc 



» 1; 



(4.68) 
(4.69) 

(4.70) 



This solution satisfies the normalization condition: Z ( cxd, u = !) = !. One can see that the main role that 
plays the second term in Eq. ( 4.67 ) is to regularize the singularity at u ^ 1 which occurs if we neglect this 
term. It should be stressed that the vertex for 2 — > 1 transition is positive for n — > 1. Eq. ( 4.70| ) suggests 
that the most important region of u for the solution is n — > 1 where Eq. ([4.67|) has a simple form 



dz_ 
dy 



u{l — u) I Kl 



dz_ 

du 



Liu) 



(5n)2 



(4.71) 



where y = r(l ^ 2) y and 

L{u) = u{l+u) + K2 (1 — u) = < 

First, we go to Mellin transform 



2 for (1 - n) < 2/k2 ; 
K2 (1 — n) for (1 — n) > 2/k2; 



(4.72) 



and equation for Z{oj,u) has the form 



a+ioo 



duj 



e^^ Z{uj,u) 



CO 



-^u + — L{u) 



u{l — u) 

where, searching the semi-classical solution, we assume that Z{uj,u) = exp {^(uj,u)). 
Eq. ( 4.74| ) has two roots: 



(4.73) 



(4.74) 



± 



2 L{u) 



1 ± Jl + 



4 uj L(u) 
Ki u{l — u) 



(4.75) 



The initial condition is Z{y = 0, n) = u. The general solution to Eq. ( 4.74 ) has the following form 
Z{y,u) == ^ e-^|e*^(-'")0+(.;) + e*-(-'")<A-(^)| 

Ja-ioo 27ri I J 



(4.76) 



where functions and should be found from the initial condition at 3^ = 0. 

The asymptotic behaviour of the solution is determined by the region of small to . One can see from 



Eq. ( [4 .761 ) that the first term in this equation at small lo has the form of 



z^+\y,u) 



^<S>{y=oa,u) 



a+ioo 



27ri 



(4.77) 



and has a chance top approach the asymptotic solution discussed above. The second term for cj <C 1 leads 
to 



z(-\y,u) = r ^ e-y le-^^'<^H r(^) 

Ja-ioo 27rz L J 



and its behaviour depends on the type of singularity at w = in function (co). For 
solution generates a constant contribution. We will try to find the solution considering 



Therefore, we want to find 



from the initial condition 



Z(+)(3^ = 0,u) = e*(^=°^'") / 

J a 



27ri 



r(^) 



(4.78) 

'{uj) = I/lj this 
(a;) = 0. 



(4.79) 



The main observation is that ^^{u!,u) (see Eq. ( [4.75D ) in Eq. ( [4.79| ) can be taken in the form 

Ki I / 4ki lo Liu') 1 



du 



1 + 



u' 1 - u' 



(4.80) 



smce 



2 L(n') 

'^^"^ » 1 either because u is close to unity or a; « — » 1. Using Eq. ( 4.72 ) we can take the 



Kl U (1 — u) 



integral over u' in Eq. ( 4.80 ) and reduce to the form 

2 



— (1 - n G 

2 \K2 



1-n 



(4.81) 



+ 



2k2 



ln(l - n) + 



K2 



In 



1- 



e 



:i-n)-- 



It easy to see that we can choose to satisfy the initial condition of Eq. ( 4.79 ) the function (p^ico) in 
the following form 



^-'(u;) = - + 



ijj 



(4.82) 



+ 



r(f 



laio] r (§ + J a uS) / a , 



r(a) 



r (a — \Ja lS) r (a + y^cFoJ) 
r(2a) 



where a = k,i/k,2- 

Indeed, for [1 — u) < 2/k2 only the first term in (/'^(cj) contribute leading to Z{y = 0,u) = 1. This 
result is the same as Z{y = 0,u) = u within our accuracy since u = 1 + 0(1/ K2)- 

The rest of Eq. (|!8|) gives a small contribution of the order of e "2 . For (1 — u) > 2/k2 the 
contribution of the first term cancels the contribution of the last one, while the second term leads to 
Z{y = 0, u) = u. Each T function has a series of poles at ^auj = a/2 + / ( / is integer number ) with 
the residues (—1)'//!. The sum over I for the first term in 2F\ looks as 



(i-n)-t5: 



-1)^ r(a + /) a: , 



i=0 



r(a) 



a T 2 



(4.83) 



with r 



. Taking into account all three terms in 2-^1 we have 



Z(+)(3^ = 0,u) = (l-u)-t 



a T 2 



(1 + r)" (l + r)2 



n 



(4.84) 



The behaviour at large y in Eq. ( 4.76| ) comes from the region of small lo. The first term stems from 
the pole in uj which leads to the asymptotic solution. The energy corrections comes from the singularities 
at negative a; in F functions. Using the following formula (see 8.326(2) in Ref. 



F(f 



■) F (f + 



F(a)F(f) 



fc=0 ^2 



(f + A:)2 + aw 



(4.85) 



we can see that the solution steeply approaches the asymptotic one having corrections that are proportional 



to e 2 



y 



A C\ 



Eq. ( 4.82| ) shows nicely the general property of the Sturm-Liouville equation [49|: the discrete spectrum 
at negative uj. It should be mentioned that our equation can be reduced to Sturm-Liouville equation in a 
general case. We also would like to stress that the semi-classical approach is valid since from Eq. ( [4.75 ) 



one can see that 1 ^ <I>t, 



Concluding this discussion we would like to summarize that we found the analytical solution in the 
entire phase space of our variable Y and u. 



Coming back to discussion of Pomeron interaction, we see that Eq. ( |4.67| ) leads to a different interac- 
tions of the Pomerons if we write the equation for the functional N of Eq. ( |3.32| ) . The number of Pomeron 
vertices that will be generated by Eq. ( [4.67| ) is greater than in the previous version of the generating 
functional Z suggested in Ref. . Indeed, the equation for N looks as 



dN 

w 



f)7 rfl N 

r(l-2)7(l-7)^ + {r(2-l)7(l-7)'(2-7) +r(2-3)7'(l-7)}^ 



One can see that Eq. (4.86) the following vertices for Pomeron interactions: 



(Pomeron intercept) P 
P - 
2P - 
2P - 
2P - 
2P - 



^ P = r(i ^ 2) 

2P = -r(1^2) 

> P = 2T{2 1) 

2P = r(2 ^ 3) - 5r(2 ^ 1) f 

3P = - r(2 ^3) + AT{2^ 1) 

4P = -r(2^1) 



r(2 ^ 3) 

i r(2^3) 



(4.86) 



(4.87) 
(4.88) 
(4.89) 
(4.90) 
(4.91) 
(4.92) 



One can see that r(2 — > 1) leads only to small corrections for 2P — > 2P and for 2P 3P transitions 
but generates correctly the 2P — > P vertex. 



It should be stressed that Eq. ([4.67] ) leads to quite different noise term in Eq. (3.10) which looks as 
follows 



dY 



with the noise term 



r(i ^ 2) ($(y) - ^\Y)) + c{Y) 

<\aY)\>= 0; 



(4.93) 
(4.94) 



< \C{Y) C{Y')\ >= 2 (r(2 ^ l)$2(y) |i _ $2(y)| j_ r(2 ^ 3) <^{Y) {1 - $(y)f j 5{Y - Y') . 
The noise term in Eq. ( [4.94[ ) is determined by the 2^1 transition only at 1 — ^(Y) ~ al <C 1, while in 



the entire kinematic region outside this band the dynamics is due to the second term in Eq. (4.94). 



4.5 A practical way to find solution: Monte Carlo simulation. 



In this section we consider the BFKL Pomeron calculus in the form of Eq. ( p.l8 ) which leads to the simplest 
approach in the framework of the generating functional aopproach. This approach is not onlapproachplest 
but also it is free from all troubles related to the negative contribution for the process of 2P — > P transition. 
Repeating procedure discussed in section 3.2 one can see that Eq. ( 3.18 ) for action leads to the following 
equation for the generating functional Z{Y, [u{b, k)]) defined as 



(4.95) 



/n 
Pn{Y - Yo;bi,ki;...; hi, h; . . . ; bn, kn) JJ u{bi, h) S Xi <f k 

n=L i=l 



where u{bi,ki) are arbitrary functions. The equation has the form 



dZ {Y -Yo;[u{b,k)]) 



dY 



X[u{b,k)] Z {Y-Yo;[u{b,k)]) 



with 



X[u\ 



I 



dHd'k 



2tt 

4:TT''a:s\^ as 



k' K{k, k')u{b, k') ^ + n(6, k) u{b, k) 



Nr 



6u{b, k' 
1 6^ 



6u{b, k) 



2^Mb,k)u{b,k)-u{b,k)) 



(4.96) 

(4.97) 
(4.98) 



Eq. (4.97) and Eq. ( [4.98 ) show that the 2P — > IP transition can be written as the two dipole to 
one dipole merging with a positive probability. Therefore, Eq. ( 4.96| ) has a very simple probabilistic 
interpretation which can be written as the following Markov's chain: 



d Pn{Y] ■ ■ ■ ]bi — b, ki — k] . . . bn — b, k^ — k) 
dY 



(4.99) 



+ 



^ ^ [Pn—i{Y; . . . ;bi — b,ki — k] . . . ; bn—i, k^—i] 



d k K{k, k ) Pn{Y', . . . ]hi — b,ki — ; . . . ; Xn — 6, kn — fe) 
4 7r Qs\ as 



27r 



{Pn+i{Y; . . .■,bj = b,kj = k;. . .;bi = b,ki = k; . . .■,bn = b,kn = k)- 



i>j 



^ni^'i • ■ ■ ]bj — b, kj — k] . . .bi — b,ki — k; . . . ;bn — b,kn — )) i 



(4.100) 



(4.101) 



This set of equations can be solved numerically and it gives a practical way to discuss the influence of 
the Pomeron loops on the solution for the scattering amplitude at high energies. 



A O 



5. High energy amplitude in QCD: solution to the functional equation. 



In this section we will find the solution to the general equation of Eq. ( p.29 ) using the strategy suggested 



in Ref. [^]. This strategy directly follows from the solution for the toy model and consists of three steps: 
1. We find an asymptotic solution to the functional equation |3.29 with zero l.h.s. using the semi- 



classical approach. In other words looking for solution in the form Z (Y = oo, [uj]) = exp {$ (Y = oo, [uj])} 
assuming that 

OUi OUj^ OUk OUi 

In searching this solution we can neglect the second tern in Eq. ( |3.2S| ) but take into account the decay of 
two dipoles into three dipoles^. 

2. Search for the function u{x,y) = UQ{x,y) which satisfies the equation: 

^ / nLi d^Xid^Vi (f z 2j^^ {xi,yi\z) K{x2,y2;z) {jo{x2,z) + 70(^,^2) - lo{x2,y2)) = 

= I I\i=i d'^XiSyi cP z K{x2,yi;z) 70(3:2, yi) (7o(a;i,z) + 7o(^;,y2) - 70(2:2, yi)) (5.1) 



where we assume that ■jo{x, y) = 1 — Uo{x,y) oc <C 1. Eq. (5.1) is the condition that the term in 
Eq. ( p. 291 ) that is responsible for 2 ^ 3 transition of dipoles, is equal to the second term in Eq. ( p. 29] ) 
which describes the decay of two dipoles into four dipoles. 

3. Specify the asymptotic solution Z {Y = 00, [ui]) from the normalization condition 

Z{Y = (^,[ui = uo]) = 1 (5.2) 

4. Search the solution in the form $ (Y, [ui]) = ^ {Y = 00, [ui]) + A$ {Y, [ui]) and assuming that 

f 5A<S>{Y,[u,])^'^ 



\ Su J small and can be neglected. In doing so we obtain a linear equation which we need to 

solve. 

5.1 Asymptotic solution 

The asymptotic solution to Eq. ( |3.29| ) with the term describing the transition of two dipoles to three dipoles 



has been found in Ref. |21] and we discuss it here for the sake of completeness of our presentation. This 
solution has the form 

$(y = 00, [n.]) = - a In ( J , y e{x,y) ) ^''^^ 



where a = ^ with ki and K2 are defined in Eq. ( 4.69| ) and is a function which is equal to 1 for \x\ < R 



and \y\ < R and = for |x| > i? and \y\ > R where R is the largest scale in the problem. From Eq. (|5^ 
one can see that is defined as 

J='{oo,[ji]) = I d^ xd^y e{x,y) j{x,y) and iV = / d^ x d^ y Q{x,y) (5.4) 



This term is not written in Eq. (3.29) but one can find it in Ref. |20l G 



A O 



Substituting Z{Y, [ui]) = exp$(y, [ui]) with ^{Y = oo, [ui]) given by Eq. (|57^) we have 
Kia J (f X y z K{x, y; z) ( 7(x, y) — 7(x, z) 7(2, y)) ~ ^ = ki a y xi S y\ S x\ S y2 z 

K{xi,y2;z) {-f{xi,y2)-f{x2,yi) - 7(^1, z) 7(2, j/2) 7(a^2, j/i)) Q(^^y2) ^^^^^ 

CanceUng common factors, multiplying by both sides of Eq. ( ^I^ ) and substituting = 
J x' y' 7(x', y') we obtain the following equation 

J d^ X d^y d^ X S y d^ z K{x, y; z) -i{x ,y') ( 7(x, y) - 7(3;, z) 7(z, y)) = (5.6) 

d^xi d^yi d^ xi d^ y2 d^ z K{xi,y2;z) (7(2;!, ^2) 7(2:2, yi) - 7(^1, z) 7(2, 2/2) 7(x2, yi)) 
One can see that if we denote in the r.h.s. of this equation xi x, y2 ^ y, X2 x' and y —>■ y', Eq. (|5 



becomes the identity. This means that Eq. ( ^.3| ) is the asymptotic solution of Eq. (|3.2S| ). 



5.2 7o(a;,y) 



In this subsection we will find 7o(x,y) as the solution to Eq. ( |5.lD in more direct way than in Ref. |21] 
using the explicit form of the vertex for 2 — > 4 dipole decay (see Eq. ( [4.58D ). If we change the notation in 
the r.h.s. of Eq. ( ^.ID , namely xi X2 and yi y2 one can see that 

7o(x,y) = — 7^^(x,y|z) (5.7) 

^2 

is the solution to this equation. It looks strange that 70 (x, y) has three arguments instead of two. However, 
the explicit form of 7^'^(x,y|z) in Eq. ( [4.56 ) shows that this function is actually a function of only two 



variable x — z and z — y. We can change all integration variable shifting them by z without inducing any 
change in the integral. It means that 7o(x,y) is equal 

1 x^ 

70(2^, y) = ^i"^'^ (5-8) 



Eq. ( |5.8D gives the same 7o(x,y) as was derived in Ref. 1 21] from rather lengthy discussions 



Using Eq. ( p.8D we can build the solution that satisfies the modified initial condition of Eq. (|5 
(Z(y, [yi = 7o,i]) = !)• This solution has the form 

Z{Y,hi\) = — — (5.9) 

V J^(oo; 7oi(x,y) ) 



A A 



5.3 Approaching the asymptotic solution 

Substituting <I> {Y; [7i(x, y)]) = ^ {Y = oo; [7i(x, y)]) + A<I> {Y; [yi{x, y)]) we obtain a linear equation for 
A<1> if we assume that 

SA<S> (y;[7.(x,j/)]) ^ 6'^> = oo; [7(x, y)]) , 
6j{x,y) Sji{x,y) 
and 

(5^a>(y = oo;[7.(x,j/)]) d-^ (Y = oo; [7(x, y)]) <5A$ (F; [7^(x, j/)]) 

Sl{xi,yi)6-f{x2,y2) hi{x,y) h{x,y) 

We can use the solution for the toy model (see section 4.4) to show that these assumptions are rea- 
sonable but we have to check Eq. ( 5.10| ) after finding the solution for A$. It should be mentioned that the 
second equation is valid since = for the asymptotic solution. 

The equation for A<I>(y; [7^]) has the form pi] ]: 
9A$(y;[7,]) _ f ,2^,2..,2. z^.. .,^ ,^ .a ^ ^ ^^(^5 N) 



(fx(fy(fz K{x,y;z) {'y{x,y) - 7(x, z) 7(2, y) ) 
dy J oj(x,y) 



+ 



1 f 

— / (i^Xid^yj d^z 27^^ (xi,yi|z) K{x2,y2;z) {j{x2, z) + j{z,y2) - -f{x2,y2)) 
'^^ •' 1=1 

+ 2— — - — - / S z(f xi(f yi(f X2(f' y2 K{xi,y2;z) 
J^((X),[7oi]) J 



6^'^iY = oo;[^,ix,y)] ) , ^ 1 f^^ 
6-f{xi,yi)6j{x2,y2) 



X {l{xi,y2)l{x2,yi) - j{xi,z)-f{z,y2)j{x2,yi) ) ^^^1^^' (5.11) 

oj{xi,yi) 

Eq. ( p. Ill ) is the Liouville-type equation with the source term which is given by the two dipoles to four 
dipoles decay. It is easy to reduce this equation to the Liouville equation without the source term and even 
linearize it if we assume that A7(3^; x, y) = 7(3^; x, y) — 7o(3^; x, y) <C 70(3^; x, y) where 70 is the solution 
to Eq. ( |5.1| ) and using the explicit form of $(y = 00, [7^]). The equation has the form 



dy 

(fzK(x,y;z) (A7(3^;x,z) + Ai{y;z,y) - A7(3^;x,y)) 



d^{x,y) 

which can be solved assuming that A<I>(y; [A7j]) = A<I>([A7j(y; Xj, y^)]) (see Refs. H^, ^] for details). 
For function A7(y; x, y) we obtain the equation 



d A7(3^;x,y) 



j (fzKix,y;z) {A^iy;x,z) + Aj{y;z,y) - Aj{y;x,y)) (5.13) 



dy 

This equation is just the BFKL equation but with the negative sign in front of l.h.s. 

The initial condition for A^([y(Y;xi,yi)]) stems from the equation Z (Y = 0, [ui]), that leads to 

exp{$(y = oo,[7i]) + A$(y = 0,[7i])} = l-7(y = 0,x,y) (5.14) 



A ir 



or 



. , X \ x(P-y<$){x,y) = ^\x,y) 

JF(3; = cx); [7o,iJ) 



The solution to Eq. ( ^.13 ) can be easily found using the Green function of the BFKL equation (see 
Eq. (PD), namely, 



A7(3^;x,y) = j d^ x' d^ y' G {x,y-Y\x' ,y' -Y = ^) ^^{x' ,y') (5.16) 

where A7(x',y') is an arbitrary function and (5 (x, y; y'; F = 0) is the Green function of Eq. ( ^.131) 
which satisfies the initial condition (see Eq. (120) of Ref. ||2^ ) 



G{x,y;Y\x' ,y';Y) = 6^^\x - x ) 6^^'> {y - y ) + 6^'^\x - y) 6^'^\y - x) (5.17) 
Such a Green function is equal to 

G{x,y;Y\x',y';Y = 0) = e{Y) x (5.18) 



X f; / ^d^xoe-^"--)^ . ""'.^j'^^ ,.^ Eix,y;xoW)E*ix',y';xoW) 
^ J [x — x')^ [y — y')^ 



n=— oo 



where all ingredients of Eq. ( 5.18 ) are defined by Eq. (^) - Eq. (|2.10| ). One can see that A7(3^;2;,y) 



decreases as exp {—uj{n = 0, z/ = 0) y) at large values of rapidity Y. 

Using Eq. (|1^) and Eq. (|5.14|) we find the functional A^>(y; [^7^]) 

A*(y;[A„]) = -A-,(Y;..y) - , Ll^pL^jp^ ,5.13, 

J'iY = oo; [joA) 



with A7 given by Eq. (5.16) while the final answer for the generating functional Z has the form 

Z(y,[A7,]) = exp{$(y==oo;[7o(x,y) + A7(y;x,y)]) + Acl>(y; [A7(y; x, y)])} (5.20) 

Using Eq. ( |5.20| ) we can easily calculate the behaviour of the scattering amplitude at high energies as 



it has been suggested in Ref. |l21 |. 



6. High energy scattering amplitude: semi-classical approach 

6.1 General equations for scattering amplitude in semi-classical approach. 

The method suggested in the previous sections leads to the exact solution at the ultra high energies, but it 
does not look very practical to reconstruct the behaviour of the scattering amplitude in the entire kinematic 
region starting from low energies. It is especially clear in comparison with the remarkable progress in the 
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searching of the solution for the mean field approach (Balitsky-Kovchegov equation). For the Balitsky- 



Kovchegov equation we have developed both analytical |5^, 51] and numerically^ approaches which give us 
the effective way for calculation of the scattering amplitude in the accessible kinematic region of energies. 
We firmly believe that the method of the generating functional will lead to a practical application in spirit 



of Ref. 1 53 1, namely, to creation of the new Monte Carlo codes which will allow to solve the equations 



for the scattering amplitude as well as to give some predictions to the inclusive observable expanding the 
region of application of this approach. 

In this section we develop the semi-classical approach to the amplitude determined by the BFKL 
Pomeron Calculus in the functional integral form of Eq. (l2.15| ) or Eq. ( |3.2D . It is well known that the 
semi-classical approach is related to the fields ^ and that satisfy the equation 

6S = (6.1) 



Making variation of S given by Eq. ( 2.16| ) - Eq. (|2.18|) with respect to both fields we obtain (see also Ref. 

i) 

dNi.,y;Y'l ^ ^^^^^ 



dY' 

as_ 
2tt 



(fzK{x,y\z) {N{x,z;Y') + N{z,y;Y') - N{x,y;Y') - N{z,y-Y') N{z,y-Y')) 



^ t Go{x,y-y\x',y'-y)^^^^K{x\y'\z) d'zN+{z,y',Y-Y')N[x',z-X) 



(27r)* TT J ^ ' {x— y') 

dN+{x,y;Y - Y' 



dY^ 



(6.3) 



1^ j cfzK{x,y\z) {N+{x,z;Y -Y') + N+{z,y;Y -Y') - N'^{x,y;Y -Y')- 
N+{z,y■Y -Y')N+{z,y;Y -Y')) 

^ Go{x,y;YV,y'-,Y') ^-^^,K{x',y'\z) d'zN{z,y'-,Y') N+{x',z-,Y-Y') 



(27r)4 TT J ' ' ^ {x'-y'Y 

where we use Eq. (|2^ ) and the new normalization for $ and <I>^, namely, <^(x, y; y) = N{x,y;Y)/4:TTas 
and ^~^{x, y; Y) = N'^{x, y; Y)/4TTas, which corresponds to the normalization of the scattering amplitude. 



The mean field approximation follows from Eq. (6.2) and Eq. ( |6.3| ) if we neglect the last terms in these 



equations which describes the two Pomeron to one Pomeron merging. The boundary conditions for these 
equations are (see Eq. ( |2.18| )) 

N{x,y;Y' = Yo) = Anas Ttar{x,y); and N+{x,y;Y - Y' = Yq) = An as Tpr{x,y) (6.4) 

Prom Eq. (|6.4D one can see when the mean field approximation could work. It should be either large T^^r 
as in the case of scattering off nuclei for which Ttar cx; Aa and N{x,y;Y' = Yq) oc asAa ps 1[^]; or the 
size of the target dipole is much large than the size of the projectile and , therefore as^tar ^ cts,pr 0]- In 
these both cases the interaction of the initial field N{x,y;Y' = Yq) leads to large contribution and the last 
term can be neglected. 
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6.2 Linear approximation to the general equation and Mueller-Shoshi band. 



As it was shown I], H, H, H that we can find the typical scale of the process (so called saturation 
momentum Qs without finding the explicit solutions of non-linear equations (see Eq. ( |6.2D and Eq. ( |6.3[) ). 
The equation for Qg is actually a value of {x — y)^ « l/Q^ at which 

N^^^\r' = {x-yf = r', = l/Q',;Y') ^ 1, (6.5) 



where N^^^^ is the solution to the linear equation (see Eq. ( |6.2| ) and Eq. ( |6.3| ) but without non-linear 
contributions). 

This equation has been solved (see Refs. [|l|, |5^) and the following expressions for the saturation 
momentum has been derived 



Q^(r ; fixed as) = Ql{Y,) exp ( {Y' - Yo) ) 



(6.6) 



Q2(y;running as) = Ql{Y,) exp ( ^j'll^^^^^ - T^^) j ; (6.7) 



where the value of 7cr can be found from the equation |jl], [5^ : 



(6.5 



1 - 7cr d'-fcr 

and 

u^il) = asxi-f) = as (2V(1) - i^il) - Hi - l) ) (6.9) 

In the mean field approach based on the Balitsky-Kovchegov equation for dipoles which size are less 
than the 1/Qs we can safely use the linear BFKL equation to obtain the scattering amplitude. As one can 
see from Eq. ( |6.2D and Eq. ( |6.3D that we need to demand that the amplitude N~^(x,y;Y — Y') should be 
smaller than unity {N~^ (x,y;Y — Y') < 1). This condition means that the dipole size x — y should be 
smaller that a new saturation scale Qs(Y — Y') which we can obtain from the same equations Eq. ( pi 
and Eq. (|3) but changing Yq ^ Y and (y' = Yq) Q^^{Y -Y' = Yq). 



In Fig. 23 we plot the kinematic region ('band') in which we can replace the system of non-linear 
equations of Eq. ( |6.2| ) and Eq. ( |6.3| ) by the BFKL linear equation. Two lines in this picture show the 
two saturation momenta: Qs{Y') ans Qs(Y — Y') in the case of fixed QCD coupling, for the amplitude 
of two dipole scattering. One dipole has the value of rapidity Y' = Yq and the size xq — yo = Rq while 
the second dipole has rapidity Y and its size is equal to i? <C Rq- Therefore, for Qs{Y') the value of 
Qs{Y' = Yq) ~ I/Rq and for all sizes of dipoles smaller than l/Qs{Y') we can neglect the non-linear 
contribution that are proportional to N'^{r,Y'). This condition can be written as follows 

HR'o/r') > P^{Y'-Yo) (6.10) 
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Figure 23: The kinematic region, where we can replace the non-linear system of equations given by Eq. (6.2) and 
Eq. ( |6.3|) (see also Eq. ( ^.13| ) and Eq. ( |6.14 )), by the linear BFKL equation. i? and Rq are the sizes of projectile and 
target dipoles, R <^ Rq). 



These sizes are to the right of the left line in Fig. ^ For Qs{Y — Y') we have a different normalization, 
namely Qq{Y -Y' = Yo) = l/R^ and instead of Eq. (l6lc| ) we have 

ln(rVi?') > ^^^^ {Y -Y' -Y^) or HRl/r^) < HRI/R") - (y -Y' -Yo) (6.11) 

1 - 7cr 1 - 7cr 

Finally, the BFKL equation can describe the scattering amplitude in the band shown in Fig. ^ which 
has the form 

\n{RlQl{Y -Y')) > HRl/r^) > ln{RlQl{Y')) where (6.12) 
In {RlQ'AY - Y')) = HRl/R') - (Y-r- Yo) ; (6.13) 

ln(i?gQ?(n) = ^^r->o); (6.14) 



This result has been derived by Mueller and Shoshi (see Ref. ||2J]) using quite a different approach 
based on the completeness relation for the dipole number density (or in other words, a unitarity constraints 
in i-channel) . In this approach the same conclusions stem from the semi-classical equation for the scattering 
amplitude. 



The most interesting result of Ref. [24| is the appearance of a new saturation scale which has the form 
^r^(RlQl{Y)) = \n{RlQl{Y)) + f =asC.Y (6.15) 

where Cg is a constant which is determined by the exact value of the amplitude where the saturation is 
declared to occur. Constant Cg is equal to 

Cg ^ 2^ + (6.16) 
I -lev 8 1n3(l/as) 
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6.3 Solution in the saturation region 



In the system of Eq. ( |6.2| ) and Eq. ( |6.3|) we have used the relation of Eq. (|2.26| ). However, we need to 
return to the general form of the action (see Eq. ( |2.17| )) and to a general form of Eq. ( |6.2[ ) and Eq. (|6.3| ). 
This general form can be written as follows 



dN{x,y-Y') 



dY' 

as 
1 as 



(6.17) 



(fzK{x,y\z) {N{x,z-Y') + N{z,y;Y') - N{x,y;Y') - N{z,y;Y') N{z,y;Y')) 



(27r)4 vr 
dN+{x,y;Y - Y' 



Go{x,y;Y'\x',y';Y') K {x' ,y'\z) <f z (L,y iV+ (z, y', y - y')) N{x',z-Y') 



dY' 



(6.18) 



as 
2tt 



(fzK{x,y\z) {N+{x,z;Y -Y') + N+ {z,y;Y - Y') - {x,y]Y - Y') 

N+{z,y;Y -Y')N+{z,y;Y -Y')) 
1 as 



(27r)4 TT 



Go{x,y;Y'\x',y';Y') K {x' ,y'\z) z {L,,y,N {z,y';Y')) N+ {x' , z;Y - Y') 



One can see that in the kinematic region where (see Fig. 24) 

.2 



and 



^Ql{Y-Y') < 1 



(6.19) 







both amplitudes N (r, Rq;Y') and N {R,r;Y — Y') are deeply in the saturation region. We can find the 
solution of Eq. ( |6.17| ) and Eq. ( |6.18| ) in this region replacing both by = 1 + AA^ and noticing that 
constant N'^ in Eq. ( |6.17| ), as well as constant N in Eq. ( |6.18 ), does not give a contribution we see that 
the asymptotic solution is = 1 and for AA^ we have the following equation 



9(L_AAr(x,y;y')) 



dY' 



(6.20) 



as 
'27r 



d^zK{x,y\z) {L^,yAN{x,y;Y')) 



as 
2tt 



d^zK{x,y\z) {L,^,AN{x,z;Y')) 



In Eq. ( |6.20D we neglect all contributions of the order (AA^)^, AA^ AA^+ and {AN+f. 



Eq. ( |6.20 ) has a solution if we assume that AA^ has a geometrical scaling [gOj, namely, AA^ is a function 
of one variable 



C = ln(^-^Q^(y')J = e - ln(i?^Qf(y' 

where Qs is given by Eq. ( |6.15 ). 

Using Eq. ( |6.6[ ) we can rewrite Eq. ( |6.20| ) in the form 



(6.21) 



x,y;z) AN{x,z;Y') 



(6.22) 




K = ln(R|/r2 



Figure 24: The kinematic region where both amphtudes N{r'^, Y') and N{r'^, E?\ Y — Y') are in the saturation 
region. 



where iV = Lx^yN{x,y\Y). In Ref. [^] it is shown that the main contribution in the saturation region 
stems from logarithmic integration, namely, 



Jl/QJY') TT 



l/Qs{Y') 

Eq. (|]2|) reduces Eq. (|a2^ ) to 



K{x,y;z) 



i/Qs{Y') (x-z)^ 



+ 



dz' 



(6.23) 



AiV(C) 

dC 



2 CAiV(C) - / dC'ANiC' 



(6.24) 



Introducing a new function AAA(C) = (iC'AiV(C') we can rewrite Eq. ( |6.24D in the form 

CsAM\,^ = -2CAM\{C) - AAA 
the solution of this equation is equal to (see 9.255(2) in Ref. |48] 



(6.25) 



AM{0 = exp 



2 a 



D 



(6.26) 



where Dp{t) is the parabolic cylinder function (see 9.240 - 9.255 in Ref. [^S]). Eq. ( |6.26| ) leads to 



AAA(C) 



exp 



(6.27) 



In semi-classical approach Eq. ( 6.27] ) leads to AA^ oc exp i— IjiJ""^-^ i^i the case if we neglect the 



difference between Qs and Qg (see Eq. ( 6.15| ). This result is very surprising. Indeed, the approach to 
the asymptotic solution in the case of the Balitsky-Kovchegov equation gives AA^ oc exp ( — 
0[ which is slower than Eq. ( |6.27| ). The estimates based on the lancu - Mueller approach |l^ gave 



r -1 



even smaller coefficient in the exponent than for the Balitsky-Kovchegov equation. The solution of the 
functional equation (see previous section and Ref. |^] ) also gives a coefficient which is in two times smaller 
than in the lancu - Mueller approach. In the last two approaches the t-channel unitarity has been heavily 
used. Certainly, the BFKL Pomeron Calculus satisfies the t-channel unitarity constraints, but we cannot 
guarantee that the semi-classical approach respects t - channel unitarity. 

The second problem with the solution of Eq. ( |6.27| ) could be our assumption that the solution in the 
saturation region has the same geometrical scaling behaviour as the solution to the Balitsky-Kovchegov 



equation. In Ref. |24| it is shown that the geometrical scaling is not correct assumption in the vicinity of 
the saturation scale in the Mueller-Shoshi band. What happens inside the saturation region we have not 
studied yet. To find out whether the geometrical scaling behaviour is an inherent property of the system 
we need to solve Eq. ( |6.2[ ) and Eq. ( |6.3[ ) in the vicinity of two saturation scales taking into account the 
non-linear corrections. We will do this in further publications. 



7. Conclusions 



We show in this paper that the BFKL Pomeron Calculus in the kinematic region given by Eq. ( |1.3| ) has 
two equivalent descriptions: (i) one is the generating functional which gives a clear probabilistic inter- 
pretation of the processes of high energy scattering and provides also a Hamiltotian-like description of 
the system of interacting dipoles; (ii) the second is the Langevin equation for directed percolation (see 
Eq. ( 3.2C| )). We are able to prove this equation for the dipoles in the momentum represetation if the 



impact parametrepresentationlarger than b ^ 1/k where k is the transverse momentum of a dipole. 

In other words, the BFKL Pomeron Calculus can be considered as an alternative description of the 
statistical system of dipoles with different kinds of interactions between them. 

One of the results of this paper is the understanding the two Pomerons to one Pomeron merging in the 



framework of the dipole approach. It turns out that the Feynman diagrams of Fig. Qq, Fig. 19 and Fig 



describe the process given by Eq. (|]5^. In this process two dipoles are produced in the first stage and 



two other dipoles in the second stage of the processes. All these four dipoles can interact with the target 
leading to Eq. ( 4.63| ) for the contribution to the generating functional. At first sight, Eq. ( 4.63| ) contradicts 



the probabilistic treatment due to the negative sign in front of this term in Eq. (3.29). Indeed, it seems 
strange that two dipoles decay into four dipoles with negative amplitude since this amplitude looks rather 
as probability of this decay. We showed in section 4 that this sign minus has a very simple explanation 
and stems from the fact that the decay of two dipoles into four actually goes in two stages. During the 
first one, one dipole decays into two with positive probability. During the second stage of the process the 
second dipole scatters off the produced two dipoles. This rescattering does not produce the new dipoles 
but change the two dipoles, produced during the first stage, to two new dipoles. This rescattering process 
goes in the Born approximation of perturbative QCD but it generates almost a pure imaginary scattering 
amplitude. This fact results in extra minus sign. The calculations of section 4 suggest the educated guess: 
all processes of merging for Pomerons ( say nP — > mP with m < n, for example such as 3P — > 2P) will 



lead to the term in Eq. (|3.29|) and Eq. ( 3.30 ) of the fohowmg type: 



2n ^ X 

It should be stressed that the vertex in front of Eq. can be found directly from the JIMWLK 

approach as as it was demonstrated in Ref. |4^]. Therefore, the evolution equations (Markov's chain) 
for the generating functional together with the JIMWLK approach lead to the selfconsistent and complete 
theory of high energy interaction in QCD. 

It should be stressed that the negative sign does not change the fact that the high energy amplitude 
can be written in terms of probabilities to find the given number of dipoles of definite sizes (see Eq. (3.25)). 



It was shown in section 3.3 that the equations for generating functional can be rewritten as the path 
integral for the partition function using which we can introduce the thermodynamic potential to describe 
the system of dipoles. 

It is shown that the question about negative amplitude does not arise if we treat the system of dipoles 
in the momentum representation. Markov's chain of equation for this system is written in the paper (see 
Eq. O)) and can be considered as a practical way to find a solution in accessible range of energies. 

The result of this paper is the solution to the linear functional equation for the generating functional 
in the high energy region. We found the semi-classical solution for the toy model in the entire kinematic 
region while we were able to calculate the asymptotic solution and the first correction to it at high energy 
in the general case of interacting dipoles. This solution shows the two BFKL Pomerons to one BFKL 
Pomeron merging is essential only in the limited region of the phase space where the resulting vertex is 
positive. 

The third result of this paper is the first attempt to solve the semi-classical equations for the BFKL 
Pomeron calculus. We confirm the Mueller-Shoshi band in ln(l/x) and ^ = ln(i?Q/r^) plane where the 
solution to the linear BFKL equation can describe the scattering amplitude which follows from the set of 
equation in the semi-classical approach. We found the solution deeply in the saturation region where the 
scattering amplitude approaches unity {N — > 1). It turns out that approach is steeper for our case than 
for the mean field approximation. 

Being elegant and beautiful the BFKL Pomeron Calculus has a clear disadvantage: it lacks theoretical 
ideas what kind of Pomeron interactions we should take into account and why. Of course, the Feynman 
diagrams in leading ln(l/x) approximation of perturbative QCD allow us, in principle, to calculate all 
possible Pomeron interactions but, practically, it is very hard job. Even if we calculate these vertices 
we need to understand what set of vertices we should take into account for calculation of the scattering 
amplitude. In this paper we show that the the BFKL Pomeron calculus together with the JIMWLK 
approach allows us to build the selfconsistent and complete theory of the high energy scattering in QCD. 
However, we lack the formalism how to extract the multi Pomeron vertices from the JIMWLK approach as 
well as the general approach which will allow us to take into account all possible multi Pomeron vertices. 
This is the reason why we need to develop a more general formalism. Fortunately, such a formalism has 
been built and it is known under the abbreviation JIMWLK-Balitsky approach EH, 37, 28 1. In this approach 



we are able to calculate all vertices for Pomeron interactions as it was demonstrated in Ref. Bql and it 
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solves the first part of the problem: determination of all possible Pomeron interactions. However, we need 
to understand what vertices we should take into account for calculation of the scattering amplitude. We 



hope that a further progress in going beyond of the BFKL Pomeron Calculus (see Refs. |45, 46|) will lead 
to such a development of the BFKL Pomeron Calculus that we will have a consistent theoretical approach. 
Hopefully this approach will be simpler than Lipatov's effective action [47| which is not easier to solve than 
the full QCD Lagrangian. 

It is well known that the mean field approach to our problem which includes only one dipole to two 
dipoles decay in the master equation (see Eq. (3.29)) has been studied quite well both analytically fSC 



and numerically |52]. We firmly believe that the probabilistic interpretation gives a practical method for 
creating a Monte Carlo code in spirit of the approach suggested in Ref. |^3|. This code will allow us to 
find a numerical solution to the problem and to consider the inclusive observables. This extension is very 
desirable since the most experimental data exist for these observables. We also hope that the semi-classical 
approach for searching the solution to the master equation will be useful for developing both analytical 
and numerical approaches which will result in predictions for the LHC. 
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A. Calculation of Go (iCi, £C2 1 a;^, 332) . 
The solution of the BFKL equation is given by 

G{x„X2,x„X2\u) - K + (n-l)V4][.^ + (n+l)V4] u:-ui.,^.) ^^'^^ 



where G,yfj_{xi,X2,Xi,X2) is the Mellin transform of Eq. ( |2.8D and Xi are two-dimensional vectors in 
complex coordinates 

Xi — Xi^x ~t~ iXi^y X^ — Xi^x ^Xi^y (A. 2) 

The function uj{u, fx) is the eigen value of the BFKL equation given by Eq. ( 2.1C| ). 



The four-point Green function is presented in terms of the hypergeometric functions [^, |6| 

Gp^i^Xl, X2, X^, X2 ) = Cix^x*^ F(h, h, 2h; x)F{h, h, 2h; x*) 
+C2X^-^x*^-'^F{l -h,l-h,2-2h; x)F{l -h,l-h,2-2h; x*) (A.3) 



r- A 



with h=^ + ii'+^,h = ^ + ii/ — ^, and x is the anharmonic ratio 



xuxr2^ (A.4) 

Xn'X22' 



Coefficients Ci and C2 are given by |2£] 



with 



,^4i,T{-iu+{l + \n\)/2)T{iu + \n\/2) 



r(ii/+(l + |n|)/2)r(-iz/ + |n|/2) ^ ^ 

As we have discussed, the high energy asymptotic behaviour stems from n = term in Eq. ([A.l|). The 



initial condition for Eq. (2.8) at y = Iq is given by the fohowing expression 

G'o(xi,X2;x'i,x2) = (A. 7) 

Gufj.\Xl, X2, Xi, X2) 



27TiJ_^ (1/2 + 1/4)2 io-io{u,n = 0) {v'^ + l/'^Y 

This integral can be taken by closing contour of integration over singularities of the integrand. 
The function 

,,2 



(j/2 + 1/4)2 



(A.J 



has two poles at | and — |. The four-point Green function of Eq. ( |A.3| ) consists of two terms, one with 
Ci(a;a;*) 2^**^ and the other one C2(xx*) 2"*^^. For small |x| this terms should be integrated closing contour in 
upper and lower semi-planes respectively. The resulting contour in the lower semi-plane runs anticlockwise 
and thus the value of the contour integral enters with a minus sign. 

The terms could be expanded in the vicinity of their poles. Let us consider the first term. We expand 
the function Ci in the vicinity of | 

' 2 r(-ii/ + i)r(i + w) (iz/ + i) ^ ■ ^ 

The hypergeometric function can be written as a sum 

F{a, b,c;x) = l+ \' ^ ^ , -r A.IO 

r(a)r(6)^ r(c + n) n! 

In the case of F{h, h, 2h; x) the singularity of T function at | can factorized out the sum 

T{2h) ^T{h + n)T{h + n)x'' 



^ ' ' ' ^ r(/i)r(/i) ^ r{2h + n) ni 



n=l 

00 



1 1 ln(l — x) , . X 



pr r 



At this stage the first term of the integrand of Eq. ( [A.7| ) can be written as 



(A.12) 



vr. 



T{-iu) 



uj{w + \f{iu-\f2 T{-iu+\)T{l+iv) 



{xx*) 



^ l /n(l-x*) 



\+iv\-, lln{l-x) 



It is clearly seen that the term of zero order in — x| has a third order pole at v = |; the term of 
first order in ln\\ — x\ has a second order pole at = |, and, the term of second order in ln\\ — x\ has a 
simple pole at = |. The contributions of those terms in the contour integral are found to be 



[-2 + ln{xx*){-A + ln{xx*))] 



— [-2 + 7 + ln{xx*)]ln[{l - x){l - X*)] 



ITT 
4UJ 



ln(l — x)ln{l — X*) 



(A.13) 

(A.14) 
(A.15) 



respectively. 

In a similar way we may expand the second of the integrand in vicinity of its pole at v = — |, namely 



(A.16) 



TT 



Tdiy) 



u;(i-ii/)3(ii/ + i)2 2 r(jz/ + i)r(l - izy) 



(xx*) 



lln{l-x*) 



The integration is performed on the lower semicircle and results in overall minus sign of the integral. The 
contributions corresponding to Eq. ( A.13 ), Eq. ( A.14 ) and Eq. ( A.15| ) are 



ITT' 



[-2 + ln{xx*){-4 + ln{xx*))] 



-[ 



-2 + 7 + ln{xx*)\ln{{\ - x)((V - x*)\ 



H /n(l — x)ln(\ — X*) 

respectively. 

Comparing the contributions we note that that of zero and second order 
in /n[(l — x)(l — x*)] are exactly canceled out, and we are left with 

-2— [-2 + 7 + ln{xx*)\ln{{\ - x){\ - x*)\ 



For small Ixl this can written as 



47r2 



ln\x\ /nil — x\ 



(A.17) 
(A.18) 
(A.19) 



(A.20) 



(A.21) 



or 



47r' 1 , , 
H In- — - m l — x\ 

U! \x\ 



Going back to complex vector representation of x and rewriting Eq. ( A.23| ) as 

Go{xi,X2;x[,X2) 



47r2 , 
H In 


XU/X22' 


In 


XU'X22' 






iO 


Xi2'Xy2 




X12XV2' 



(A.22) 



(A.23) 



we see that we reproduce the result of [26]. Therefore, we demonstrated that Eq. ( A.23| ) gives a correct 
initial condition for searching the scattering amplitude at high energies restricting ourselves by the one 
term in Eq. (p.Sp with n = 0. 



B. The path integral formalism for the generating functional. 

We want to develop a path integral formalism similar to that we found for the toy model, but where the 
probabilities to find n dipoles depend of dipole sizes. In our notation we denote by Latin index rapidity 
interval, and Greek indices relate to a size of dipole. 

As in Chapter 3.3 we introduce the creation and annihilation operators 



a{q) 



5u{q) 



a\q) = u{q) 



(B.l) 



with commutation relations [0(5), a^(g')] = 5{q — q') at fixed Y. The expression for the coherent states in 
this case takes form of 



(B.2) 



with 



(B.3) 



First we consider discrete dipole sizes qa = {L/N)a, where L is a maximal possible dipole size, N a 
number of intervals of L discretization, and a is an integer number running from to N. In this case the 
commutation relations become [a(qa)Ta^{q/3)] = 5a,i3- 

The unit operator can be written in terms of the coherent states 



n 



\<l){qa) >< (l)iqa)\ 



(B.4) 



r 17 



The operator H, defined in Eq. (|3.30| ) can be written as 



nla), a] = - ^ ^ Yy^i^2{qi3 ^ ^7 + 9A){-aU9/3) + a) {q^)a) {qx)]a{qp) 
/3 7 A 

- V2--.i{q^ + qx^ qp){aHqj)a^ (qx) - a^qi3)}^a{qj)a{qx)] (B.5) 
Following the logic of Chapter 3.3 we consider a matrix element 

|n < </'.+i('?"')i| (1 + WAY) |i n</.,(Qj >| (B.6) 

First, look at the second term of the Hamiltonian 



la' ) 13 -y \ la ) 

I n < I E E E[-^i-2(g/3 - 97 + QxWj+iiQ'rWj+ii^x)<P,iQp)] \ U \<P,iQa) > \ (B.7) 

la' J/37A la J 



In Eq. ( |B.7D we used the property of the coherent states given by Eq. 

In the continuous limit Sa,i3 is replaced by S{qa — q/3), and dcp* {qa)d(j){qa) by functional integration 
/ V(j)*V(j). 

Prom here we see that rest of the calculations is similar to that of Chapter 3.3 and we end up with the 
expression for a matrix element of an operator A between states of initial Yq and final rapidity Y 

< Y\A\Yo >~ J P$+P$y4(y)e'^ (B.8) 



where 



S = J (^J <^^{q)^Hq)dq + ni<^>+ + l,-<l>)^dY (B.9) 



with the Hamiltonian given by 



+ j d^qod'^qid^q2[Vi-.2{qo ^ qi + q2){-^^{qo) + ^+(91) + ^+(92) + $+(gi)$+(92)}$(go) 

- V2^i{qi + 92 ^ 9o){^+(9i)^+(92) + ^+(91) + ^+(92) - ^+(9o)}^^(9i)^(92)] (B.IO) 



pro 



C. The exact calculations in the semi-classical approach 



At this appendix we calculate exactly exponent in semi-classical solution and compare it with it approxi- 
mations at small and large values of uj, which ware used for calculation of corresponding weight functions. 

Remind that Eq. ( [4 .751 ) 



= 1 ± Jl + ' "."^^^^ . (C.l) 

2 L(n) I y Ki ^(1 - n) J ^ ^ 

Therefore exponent in semi-classical solution <I>^ will be defined by the following integrals: 



Jo 2 L{u') V Ki m' (1 - u') 



where 



Kl 



2L(n') 
2 

a/ 6k2 - 1 - 



arctan 



1 — K2 + 2n 
a/ 6k2 - 1 - k 



arctan 



1 - 

y^6K2 - 1 - 



(C.3) 
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Kl ti' (1 — u') 



du' 



A UJ ■ [u' {I + U') + K2{1 - U')] 



2 7o n'(l + u') + K2(l-n') 
± [//2(n)- 1/2(0)] 



1 + 



Kl m' (1 — U') 



(C.4) 



where II2 defined as following indefinite integral: 



I hit) = t 



(A+ + 9,UJ^n-t)) ■ (A- + Sij!^(l~t)) / ^ + i^t + K2f2 

^ V V '1 + 46^^11 ^ ^ 

t{l-t) 



EPi 



i ■ arcsin h 



4lJ K2 t 



+ EPi 



A- 1 - t / ' 1+ 

i • arcsin h 



8UJ K2 1 - ^ \ A~ 



j ■ arcsin 



r r\ 



1 



t(l-t) 



In the last relation we used the following notations: 



(C.5) 



EF{ip,m) 
EPi{n, ip, m) 



Kl 

de 



^/T^^rn~^I?J&) 



= Elliptic integral of the first kind 



de 



(1 - n • sin2(6')) -^l-m- sm'^{9) 



= Elliptic integral of the third kind 



(C.6) 

We can see that in spite the fact that can be calculated explicitly (in radicals), it is too complicated 
in order to use it for the calculations of the weight functions (/>^(w) in Eq. ( [4.771 ) ^-iid Eq. ( 4.78 ). Never- 
theless, this exact calculation very fruitful for two purposes: (I) using exact form of we can estimate 
how good our approach for small and large uj (or in another words how good we estimated spectrum of our 
solution), (II) once we know our spectral functions (j)^{uj) we can use this exact solution for the following 
calculations. It is important to note that Il2{t) looks as imaginary functions, but on fact the imaginary 
part is closely related to regularization of singularity at u = 0, thus it is essential in some since. The plots 
below illustrates how good our approach for small and large ui, which we used for calculations of <I>^. 




Figure 25: At this plot wc compare positive branch of our exact result (thick doted line) with for 0(a;) with 



approximation for small values of uo (i.e. uj > ^), which was used in Eq. (4.77) (thin solid line) 
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